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EXECUTIVE  SUMMARY 


THE  PROBLEM 

During  the  exercise  of  strategic  simulations  and  models  for  the 
estimation  of  the  total  number  of  injuries  and  fatalities  following  a  nuclear 
exchange,  a  number  of  critical  parameters  must  be  employed  whose  exact  values 
are  unknown,  yet  whose  impact  on  simulated  results  is  considerable.  One  such 
parameter  is  the  assumed  fallout  shelter  'stay  time'  of  the  civilian 
population. 

One  standard  currently  employed  in  the  strategic  simulation 
community  is  the  set  of  effective  protection  factors  developed  by  SRI  (4), 
based  on  the  following  stay  time  assumptions: 

(1)  7  days  fully  sheltered  followed  by  14  days  partially  sheltered  for 
shelters  exhibiting  adequate  plumbing/water  facilities,  and 

(2)  2  days  fully  sheltered  followed  by  3  days  of  partial  sheltering  for 
facilities  exhibiting  generally  inadequate  plumbing/water 
facilities. 

These  protection  factors  are  incorporated,  foe  example,  by  the  Studies, 
Analysis,  and  Gaming  Agency  (SAGA)  in  the  Single  Integrated  Damage  Assessment 
Capability  (SIDAC)  for  the  estimation  of  fallout  fatalities  and  injuries  in 
various  strategic  scenarios.  However,  the  shelter  stay  times  reflected  in 
these  parameters  represent  little  more  than  best  reasonable  estimates  of 
shelter  stay.  The  use  of  these  and  other  'reasonable'  assumptions  of  stay 
generate  a  wide  variation  in  the  number  of  simulated  fallout  casualties. 

As  a  result,  a  definite  need  exists  for  the  development  of  a 
methodology  to  justify  choice  of  stay  time  assumption  and  to  reduce  the 
variability  associated  with  the  human  factor  in  large  scale  simulations. 
OBJECTIVE 

The  objective  of  this  study  is  threefold: 

(1)  First,  to  develop  a  model  of  fallout  shelter  stay  time  based  on  the 
application  of  principles  of  human  behavior  to  an  empirical 
database; 


(2)  Second,  to  provide  a  guide  for,  and  examples  of,  the  use  of  this 
model  for  the  strategic  analyst  specifically  involved  in  the 
simulation  of  fallout  casualties  in  strategic  exchanges;  and 

(3)  Third,  to  reduce  the  wide  variation  in  simulated  strategic  cas¬ 
ualties  currently  displayed  as  a  result  of  the  spectrum  of  stay 
times  employed. 

APPROACH— DEVELOPMENT  OF  THE  MODEL 

The  model  of  fallout  shelter  stay  developed  in  this  study,  denoted 
the  Attrition  Rate  Model,  is  based  on  a  statistical  analysis  of  a 
quantitative  database  of  human  response  in  disaster  situations.  Behavioral 
studies  of  human  response  to  disasters  were  systematically  surveyed  for 
quantative  estimates  of  time  spent  in  shelters  and  qualitative  data  bearing 
on  the  behavioral  profile  of  shelterees.  Three  thousand  estimates  were 
derived  from  approximately  seven  hundred  studies.  Eight  general  categories 


of  data  were  identified: 

(1) 

physiological , 

(2) 

shelter  space, 

(3) 

shelter  type, 

(4) 

warning, 

(5) 

training, 

(6) 

shelter  management, 

(7) 

evacuation  posture,  and 

(8) 

communication. 

For  each  of  these  categories,  estimates  of  percent  shelteree  attrition  were 
derived  as  a  function  of  time  since  sheltering.  Results  were  refined  by 
eliminating  data  from  incidents  not  analogous  to  a  oost-attack  environment. 
Qualitative  behavior  profiles  derived  from  questionnaires,  interviews  and 
observations  were  developed  and  used  to  interpret  the  quantitative  data. 

The  resulting  model  expresses,  for  a  variety  of  shelter  environ¬ 
ments,  the  percent  of  shelterees  leaving  the  shelter  as  a  function  of  time 
since  sheltering.  It  was  found  that,  in  contrast  with  traditional 
assumptions  which  assume  entire  population  either  completely  in  or  out  cf 
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shelters  at  any  given  time,  that  the  expected  behavior  response  of  a 

sheltered  population  is  one  of  constant  'attrition'  from  the  shelter. 
UTILITY--A  GUIDE  TO  THE  MODEL'S  USE 

The  utility  of  the  Attrition  Rate  Model  of  fallout  shelter  stay 
is  based  on  two  features  of  the  model: 

(1)  First,  that  it  provides  a  justification  for  choice  of  stay  time 

assumption  in  strategic  simulations  by  virtue  of  its  basis  in 

a  quantitative  database,  and 

(2)  Second,  that  it  provides  the  potential  for  the  reduction  in 

the  wide  variation  of  simulated  fallout  casualties  associated 
with  the  range  of  stay  time  assumptions  traditionally  employed. 

In  order  to  allow  the  strategic  analyst  to  incorporate  these 
advantages  into  current  estimates  of  fallout  casualties,  and  to  derive  full 
advantage  from  the  Attrition  Rate  Model  in  strategic  simulations,  a  guide 
to  the  use  of  the  model  in  both  pencil -and-paper  and  large  scale  simul¬ 
ations  is  developed  and  presented.  The  guide  focuses  on  both  the  generic 
step-by-step  procedure  for  the  estimation  of  fallout  casualties,  and  is 
illustrated  by  two  explicit  applications:  a  study  of  casualties  in 
Tbilisi,  Soviet  Georgia,  and  an  investigation  of  casualties  associated  with 
a  SIDAC  scenario. 
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This  report  is  submitted  to  the  Defense  Nuclear  Agency  by  the  BDM 
Corporation  under  Contract  Number  DNA-Q01-78-C-0060  and  represents  work 
conducted  during  the  period  1  June  1978  through  1  January  1979. 

The  objective  of  this  study  is  to  develop  a  model  of  fallout 
shelter  stay  times  and  to  provide  a  guide  for  its  application  in  strategic 
simulations.  The  model,  denoted  the  Attrition  Rate  Model,  is  based  on  an 
empirical  database  constructed  from  a  systematic  review  of  behavioral 
response  in  American  disasters.  Traditionally,  stay  time  assumptions 
employed  in  strategic  simulations  have  been  based  on  little  or  no  quantita¬ 
tive  data,  but  have  represented  only  reasonable  or  computationally  conve¬ 
nient  assumptions.  The  use  of  an  empirical  database  in  the  Attrition  Rate 
Model  provides  a  justification  for  the  choice  of  simulation  stay  time,  and 
thus  contributes  to  reducing  the  wide  variation  in  simulation  results 
characteristically  associated  with  the  'human  element' . 

The  utility  of  the  Attrition  Rate  Model  can  be  measured  only  by 
its  usefulness  and  applicability  within  the  strategic  simulation  community. 
For  this  reason,  the  material  in  this  report  is  presented  in  a  concise 
jjF  manner,  specifically  focused  toward  providing  the  strategic  analyst  with  a 

p  guide  for  its  use  in  simple  'pencil  and  paper'  studies  and  large  scale 

•9  simulations,  such  as  the  Single  Integrated  Damage  Assessment  Capability 

(SIDAC). 

The  principal  authors  of  this  report  are  J.  H.  Kinrich,  N.  J. 
Maresea,  and  R.  A.  Levit. 
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SECTION  1 

INTRODUCTION  TO  THE  PROBLEM 


1 . 1  BACKGROUND 

One  key  element  involved  in  the  overall  determination  of  the 
effectiveness  of  U.S.  and  Soviet  strategic  nuclear  forces  is  the  use  of 
strategic  exchange  simulations  and  targeting  models.  Among  these  simula¬ 
tions  are  the  Single  Integrated  Damage  Analysis  Capability  (SIDAC),  which 
is  a  large  scale  nuclear  exchange  computer  simulation,  and  CIVIC,  COBRA, 
READY,  and  RISK  II  (1,  8).  These  tools  allow  strategists  to  investigate 
different  scenarios  characterized  by  variations  in  the  distribution,  type, 
yield,  and  accuracy  of  strategic  weapons;  alternate  mission  types  (counter¬ 
force,  counter-value,  etc.);  differing  target  vulnerabilities;  and  other 
strategic  issues.  The  resulting  number  of  fatalities  and  injuries  sus¬ 
tained  by  the  population  through  simulated  prompt,  collateral,  and  fallout 
effects  contributes  directly  to  estimates  of  a  country's  ability  to  re¬ 
cover;  a  process  of  significant  strategic  importance  to  any  nation. 
However,  during  the  exercise  of  these  simulation  tools  for  the  determin¬ 
ation  of  population  attrition,  a  number  of  critical  parameters  must  be 
employed  whose  exact  values  are  unknown,  yet  whose  impact  on  simulation 
results  is  considerable.  One  such  parameter  is  the  length  of  time  the 
civilian  population  can  be  expected  to  remain  in  fallout  shelters  following 
a  nuclear  attack. 

Currently,  the  Defense  Nuclear  Agency  and  others  employ  fallout 
protection  factors  developed  by  SRI  (4).  These  standard  PFs  are  heavily 
dependent  on  the  shelter  stay  time  assumptions  used  to  generate  them.  SRI 
uses  two  sheltering  assumptions  in  its  study;  these  have  been  called  the 
"Seven  Day"  and  "Two  Day"  assumptions.  For  example,  developed  countries 
are  assumed  to  have  shelters  with  sufficient  food,  water,  and  plumbing 
facilities  to  allow  a  seven  day  stay  in  the  shelter.  This  seven  day  period 
is  followed  by  fourteen  days  of  partial  shelter  occupancy:  2/3  of  the  day 
in  the  shelter,  and  1/3  of  the  day  outside.  Lesser  developed  countries  do 
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not  have  as  highly  developed  water  and  plumbing  systems.  Therefore, 
shelters  in  these  countries  are  assumed  to  allow  two  days  of  complete 
sheltering,  followed  by  three  days  of  partial  (2/3  in,  1/3  out)  sheltering. 

In  order  to  provide  information  on  the  extremes  of  shelter  stay, 
some  modelers  also  investigate  two  additional  possibilities:  Indefinite 
sheltering  (100%  sheltering  until  radiation  level  reaches  zero)  and  no 
sheltering. 

These  variations  in  currently  employed  stay  time  assumptions 
support  a  wide  variation  in  resulting  population  fatalities,  yet  each 
assumption  represents  little  more  than  a  "reasonable'1  estimate  or  compu¬ 
tational  simplification  of  shelter  stay.  To  date  the  user  has  no  empirical 
basis  by  which  to  choose  a  fallout  shelter  stay  time.  A  definite  need 
exists  for  the  development  of  a  methodology  to  justify  the  choice  of  stay 
time  assumptions  and  to  reduce  the  variability  associated  with  the  human 
factor  in  large  scale  simulations. 

Determination  of  fallout  shelter  stay  time  involves  the  evalu¬ 
ation  of  human  behavior  under  stressful  situations.  As  such,  a  behavioral 
science  approach  to  stay  time  estimation  will  provide  the  needed  human 
factors  element.  This  study  represents  an  initial  effort  to  address  these 
needs  in  a  quantitative  manner. 

1.2  STUDY  OBJECTIVE 

The  objective  of  this  study  is  threefold: 

(1)  First,  to  develop  a  model  of  fallout  shelter  stay  times 
based  on  the  applications  principles  of  human  behavior 
to  an  empirical  database 

(2)  Second,  to  provide  a  guide  for  the  use  of  this  model  as  a 
tool  in  casualty  investigation,  and; 

(3)  Third,  to  reduce  the  wide  variation  in  simulated  strategic 
casualties  currently  displayed  as  a  result  of  the  spectrum 
of  stay  time  assumptions  employed. 
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In  support  of  these  objectives,  a  model  of  fallout  shelter  stay 
times,  denoted  the  Attrition  Rate  Model,  is  constructed  based  on  a  compli- 
ation  and  anlaysis  of  quantitative  stay  time  estimates  from  behavioral 
studies  of  disasters. 

This  document  is  a  guide  to  the  modeling  methodology  and  appli¬ 
cations  of  the  model  itself.  Instruction  in  the  use  of  the  model  in  both 
"penci 1-and-paper"  and  large  scale  computer  simulation  studies  is  provided 
with  an  example  application  for  each  type  of  study.  These  examples  illus¬ 
trate  the  reduction  in  the  variability  of  fallout  casualties,  while  the 
model  itself  provides  a  justification  for  the  choice  of  stay  time  assump¬ 
tions. 

1.3  REPORT  ORGANIZATION 

The  remainder  of  this  report  is  composed  of  three  sections.  In 
Section  2,  the  methodology  employed  in  the  development  of  the  Attrition 
Rate  Model  is  presented  along  with  a  discussion  of  model  extensions  and 
limitations.  A  guide  for  the  use  of  the  model  as  a  tool  in  casualty 
investigation  and  examples  of  its  use  are  given  in  Section  3.  These 
examples  comprise  first;  an  exemplary  pepci 1-and-paper  study  of  casual¬ 
ties  associated  with  a  simulated  laydown  on  the  city  of  Tbilisi,  Soviet 
Georgia,  and  second;  a  computer  simulated  investigation  of  casualties  based 
on  a  SIDAC  red-on-blue  simulation.  Finally,  Chapter  IV  discusses  addi¬ 
tional  applications  and  extensions  of  the  model. 
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SECTION  2 

A  MODEL  OF  FALLOUT  SHELTER  STAY  TIMES 


2. 1  INTRODUCTION 

In  order  to  develop  a  quantitative  model  of  fallout  shelter  stay 
times  in  a  nuclear  environment,  an  empirical  database  must  be  developed 
which  reflects  as  closely  as  possible  the  expected  circumstances  following 
a  nuclear  exchange.  From  the  perspective  of  the  sheltered  population,  this 
environment  is  characterized  by  the  shelter  atmosphere  and  the  existence 
and/or  perception  of  a  real  external  threat. 

The  only  peacetime  experience  available  which  provides  quantita¬ 
tive  data  on  shelter  stay  times  in  real  threat  environments  is  that  of 
disaster  studies.  Indeed,  a  fundamental  assumption  of  the  present  study 
is  that  a  post-attack  environment  i_s  a  disaster,  and  that  human  response 
to  a  nuclear  disaster  is  an  extropolation  of  human  response  to  natural 
disasters.  Consequently,  the  model  of  fallout  shelter  stay  time  proposed 
in  this  study  represents  a  quantitative  summary  of  stay  times  based  on  an 
extensive  disaster  database  characteristic  of  the  American  population. 
Because  this  database  is  derived  from  a  spectrum  of  disaster  shelter  types, 
and  includes  behavioral  response  to  real  threats,  it  is  felt  to  closely 
reflect  the  key  behavioral  determinants  expected  in  a  post-attack  nuclear 
environment. 

The  development  of  this  model,  denoted  the  Attrition  Rate  Model, 
by  the  use  of  a  quantitative  empirical  database  provides  a  justification 
for  its  use  as  the  "stay  time  assumption"  in  strategic  simulations.  Cur¬ 
rently  employed  stay  time  assumptions  are  based  on  little  or  no  data,  but 
represent  at  best  only  reasonable  approximations,  and  at  worst,  computa¬ 
tionally  convenient  parameters.  As  such,  there  is  no  specific  criteria  by 
which  to  choose  among  assumptions,  yet  simulation  results  are  quite  sensi¬ 
tive  to  the  assumptions  employed.  Thus,  the  Attrition  Rate  Model  provides 


a  distinct  advantage  to  the  strategic  analyst  by  providing  a  more  justifiable 
choice  of  stay  time,  and  consequently  a  reduction  in  the  variation  in 
simulation  results  traditionally  associated  with  the  'human  element'. 

2.2  OVERVIEW  OF  MODEL  DEVELOPMENT 

The  approach  employed  in  developing  the  Attrition  Rate  Model  of 
fallout  shelter  stay  times  involved  the  construction  of  a  quantitative 
database  of  stay  time  estimates  extracted  from  relevant  behavioral  studies, 
and  an  analysis  of  the  data  for  application  to  a  post-attack  environment. 
The  resulting  'model'  is  cast  in  the  form  of  graphical  displays  (with  asso¬ 
ciated  analytical  representations)  illustrating  the  percent  of  the  shelter 
ed  population  expected  to  exit  the  fallout  shelters  as  a  function  of  time 
during  the  first  ten  days  following  shelter  entrance. 

Approximately  seven  hundred  studies  of  human  response  to  disaster, 
isolation  and  confinement,  and  shelter  occupancy  were  reviewed,  from  which 
about  three  thousand  estimates  of  stay  time  were  extracted.  These  data  were 
reduced  to  five  hundred  and  seventy  six  points  by  eliminating  of  data  from 
incidents  not  analogous  to  a  post-attack  environment.  Within  this  final 
database,  a  frequency  analysis  was  performed  to  identify  general  catagories 
of  data  which  characterized  reported  motivations  of  shelter  exit.  These 
included  physiological,  shelter  type,  shelter  space,  warning,  training, 
shelter  management,  evacuation  posture,  and  communication.  Data 
characteristic  of  each  of  these  parameters  were  least-squares-fitted  to 
linear,  and  in  some  cases  simple  power  law,  response  curves.  Qualitative 
behavior  profiles  derived  from  questionnaries ,  interviews,  and  observations 
were  also  developed,  and  used  to  interpret  these  quantitative  data.  Figure 
2-1  illustrates  the  overall  approach  employed. 

8y  combining  data  representing  stay  time  estimates  reflecting  the 
best  configurations  of  each  of  the  above  categories  and  data  representing 
worst  configurations,  'best  case'  and  'worst  case'  extremes  in  expected 
shelter  exit  response  were  developed.  In  this  form,  the  model  represents  a 
useful  tool  for  the  investigation  of  the  expected  variation  in  fallout 
casualties  due  to  the  range  of  expected  human  response. 
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Figure  2-1.  Development  of  the  attrition  rate  model  of  fallout 
shel ter  stay  times . 


2.3  DEVELOPMENT  OF  THE  MODEL 

2. 3. 1  The  Database 

The  source  literature  reviewed  for  the  development  of  the  Attri¬ 
tion  Rate  Model  included  approximately  seven  hundred  studies  of  disaster, 
isolation  and  confinement,  and  fallout  shelter  occupancy.  The  overwhelming 
majority  of  the  quantitative  data  subsequently  employed  in  the  model  data¬ 
base  was  based  on  disaster  studies  of  the  National  Academy  of  Sciences-- 
National  Research  Council  (Disaster  Research  Group),  the  National  Opinion 
Research  Council,  and  the  Ohio  State  Disaster  Research  Team  (2,9,10,11). 
Table  2-1  depicts  the  nature  of  the  disaster  research  reports  included  in 
the  study.  A  decision  was  made  not  to  include  results  of  fallout  shelter 
occupancy  and  isolation/confinement  studies  since  the  component  of  behavior 
reflecting  the  presence  of  a  real  external  threat,  as  would  be  character¬ 
istic  of  a  nuclear  environment,  was  absent  from  these  studies. 

From  the  disaster  studies,  quantitative  estimates  of  shelter 
stay  times  were  identified  and  added  to  the  database.  For  this  purpose, 
the  act  of  'sheltering'  was  considered  to  be  any  overt  behavior  whose 
objective  was  to  partially  or  completely  protect  the  individual  from  the 
physical  effects  of  an  incident.  Thus,  in  the  case  of  a  flood,  for 
example,  the  'shelter'  may  have  been  a  rooftop,  or  in  the  case  of  an 
earthquake,  under  a  table  in  the  basement.  Actual  stay  time  data 
represent  first  hand  reports  via  questionnaires  and  interviews  with 
study  participants  or  disaster  victims,  as  well  as  observations  reported 
by  reliable  observers  such  as  the  Red  Cross,  police,  fire  rescue,  and 
other  personnel  involved  in  disaster  recovery. 

Stay  data  extracted  from  these  studies  were  recorded  by  including 
event  type  and  description,  shelter  type,  stay  time,  reason  for  leaving, 
and  time  of  arrival  of  aid.  All  data  were  coded  and  placed  on  computer 
cards  for  subsequent  processing.  As  a  check  on  internal  consistency, 
the  source  documentation  was  divided  in  two,  and  data  were  extracted 
from  each  half  by  different  analysts.  The  resulting  response  curves 
based  on  these  two  sets  of  data  were  found  to  correlate  closely. 
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Table  2-1.  Sources  of  disaster  studies. 


BEHAVIORAL  ASPECTS  SOURCES 

OF  FALLOUT  SHELTER  STA  Y  DISASTER  RESEARCH 


DISASTER  AGENTS 
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STUDIES 
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4 

3 
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7 
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3 

2 

19 

2 

EARTHQUAKES.  ETC. 

8 

8 

1.831 

10 

EPIDEMICS  AND 
EPIDEMIC  THREATS 

5 

5 

2.487 

7 

EXPLOSIONS 

AND  FIRES 

13 

13 

678 

11 

FALSE  ALERTS 

6 

7 

2.953 

7 

FLOODS 

12 

16 

3.319 

27 

HURRICANES 

AND  TYPHOONS 

12 

9 

364 

9 

MINE  DISASTERS 

2 

3 

297 

5 

TORNADOES 

20 

31 

2,092 

34 

TOXICOLOGICAL 

SUBSTANCES 

8 

8 

227 

6 

WORLD  WAR  II 
BOMBINGS 

4 

6 

7,163 

4 

MISCELLANEOUS 

6 

3 

18 

4 

TOTALS 

103 

114 

21.624 

121 

21  21/78W 
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The  initial  three  thousand  point  database  developed  in  this* 
manner  represented  a  summary  of  quantitative  stay  time  estimates  and  event- 
character! sties  fro m  reliable  observers,  3nd  provided  a  foundation  r'or  the* 
development  of  the  Attrition  Rate  Model.  '* 


2. 3. 2  Data  Analysis 

In  order  to  obtain  a  database  most  representati ve  of  a  post-  ] 

] 

attack  environment,  a  number  of  criteria  were  imposed  by  which  to  eliminate] 
inappropriate  data  from  the  initial  database.  These  criteria  demanded  that 
the  data  represent: 

(1)  incidents  characterized  by  brief,  intense  expenditure  of  energy; 

(2)  stay  time  estimates  before  the  arrival  of  aid; 

(3)  data  collected  within  30  days  of  the  incident;  and, 

(4)  data  collected  within  impact  and  fringe  areas. 

The  application  of  these  criteria  reduced  the  dataoase  f^om  three  tncusand 
to  approximately  six  hundred  data  points. 

Within  this  final  database,  a  frequency  analysis  was  employed  to 
identify  categories  of  data  associated  with  shelter  exit.  Eight  such  cate¬ 
gories  were  identified: 

1)  physiological,  2)  shelter  space,  3)  shelter  type,  4)  warning,  5) 
training,  6)  shelter  management,  7)  evacuation  posture,  and  3)  com¬ 
munication.  The  data  were  partitioned  among  these  categories,  and  a  ^east- 
squares  linear  fit  was  performed  to  develop  response  curves  for  each  cate¬ 
gory.  The  resulting  curves  are  shown  in  Figure  2*2.  Note  that  in  some 
cases,  it  was  possible  to  fit  the  data  to  a  simple  power  law  more  accu¬ 
rately  than  to  a  straight  line. 

For  the  purpose  of  employing  these  results  in  the  investigation 
of  fallout  casualties,  two  additional  cases  were  developed:  Best  Case,  and 
Worst  Case.  For  the  Best  Case  response,  all  data  representati  ve  of  the 
best  cases  of  each  of  the  above  eight  categories  were  included.  These  data 
were  then  fit  to  a  straight  line,  as  shown  in  Figure  2-3a.  Thus,  this 
result  summarizes  the  expected  responses  reflecting  adequate  food,  water 
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Figure  2-2.  Graphs  of  the  general  attrition  rate  equations. 


20 


Propensity  for  Attrition  (%) 


c 

o 


t 

< 


BASIC  PHYSIOLOGICAL 

Food  Availability 
No  Food  Available 


o 

>• 


c 

<D 

a 

o 


Linear  y  -  a  +  bx 
90%  Confidence  Interval 
Slope  =  0.2575 
Y  Intercept  =  0.0035 
Coefficient  of  Correlation  -  0.9323 
Standard  Error  of  Estimate  -  4  4303 
T 


t — i — r  .  . 

120  144  168  192  216  240 


Time  Since  Sheltering  (Hours) 


60-, 


BASIC  PHYSIOLOGICAL 

Water 
500  Ml/Day 


Linear  y  =  a  +  bx 
90%  Confidence  Interval 
Slope  =  0.4260 
'  Q  Y  Intercept  =  4.6555 

/  Coefficient  of  Correlation  =  0.8193 
Standard  Error  of  Estimate  =  9.0935 


24  48  72  96  120  144  1^8  1^2  216  240 

Time  Since  Sheltering  (Hours) 


V 

< 


at 

c 

« 

a 

o 


23-i 

22- 

21- 

20- 

19- 

18- 

17- 

16- 

15- 

14- 

13- 

12- 

11- 

10- 

9- 

8- 

7- 

6- 

5- 

4- 


W  ARMING 

Warning 


Power  y  =  a*x**b 
b  =  0  4601 
a  =  2  0696 
R  Squared  =  0  9379 
Standard  Error  of  Estimate  =  0.1271 
X 


48 


72 


96 


“1 - 1 - 1 - 1 - 1 - 1 

120  144  168  192  216  240 


Time  Since  Sheltering  (Hours) 


WARNING 

No  Warning 


£ 

c 

o 


< 


at 

c 

® 

a 

o 


Linear  y  =  a  +  bx 
90%  Confidence  Interval 
Slope  =  0.2069 
Y  Intercept  =  8.2309 
Coefficient  of  Correlation  =  0.9725 
Standard  Error  of  Estimate  =  3.0381 

T - 1 - 1 — T — I — I 

120  144  168  192  216  240 


Time  Since  Sheltering  (Hours) 


Propensity  for  Attrition  (%)  I  Propensity  for  Attrition  (%) 


Time  Since  Sheltering  (Hours) 


44  0-i 
41  5- 

1 

39  0- 

c 

36.5- 

o 

♦3 

34  0- 

31.5- 

£ 

29.0- 

< 

k_ 

24.0- 

o 

21  5- 

> 

19.0- 

4-» 

’</> 

16.5- 

c 

0) 

14.0- 

a 

11.5- 

o 

90- 

CL 

6.5- 

40- 

COMMUNICATION 

Present 

Low 


Power  y  a*x**b 
0  6940 
a  -  1  1850 
R  Squared  --  0  9826 
Standard  Error  of  Estimate  =  0  1157 

~ I - 1 - 1 - 1 - 1 - 1 - 1 1 

72  96  120  144  168  192  216  240 

Time  Since  Sheltering  (Hours) 


f 


SHELTER  TYPE 

Formal 
Impact  Area 


Linear  y  =  a  »  bi 
90%  Confidence  Interval 
Slope  =  0.1412 
Y  Intercept  =  0.6425 
Coefficient  of  Correlation  =  0.9758 
Standard  Error  of  Estimate  =  2.0781 

.  .  “I - 1 - 1 — T~I - 1 — I - 1 — I 

0  24  48  72  96  120  144  168  192  216  240 

Time  Since  Sheltering  (Hours) 


c 

o 


C 

< 


(A 

c 

® 

a 

o 


71 

66- 
61  - 
56- 
SI  - 
46- 
41  - 
36  - 
31  - 
26- 
21  - 
16 
11  - 
6- 
1  - 


SHELTER  TYPE 

Expedient 
Impact  Area 


Linear  y  =  a  +  bx 
90%  Confidence  Interval 
Slope  0  2938 
Y  Intercept  =  10  8708 
Coefficient  of  Correlation  =  0  9192 
Standard  Error  of  Estimate  -  7  0096 


— 1 - | - | - | - 1 - 1 - 1 - ! - - - I 

24  48  72  96  120  144  168  192  216  240 

Time  Since  Sheltering  (Hours) 


Figure  2-2.  Graphs  of  the  general  attrition  rate  equations  (continued) 
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Figure  2-2.  Graphs  of  the  general  attrition  rate  equations  (continued) 
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Figure  2-3.  Graphs  of  the  attrition  rate  equations  -  Best  and  Worst  cases 
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and  shelter  space,  and  the  presence  of  warned,  trained,  and  well  managed 
shelterees  with  adequate  communication.  The  Worst  Case  is  depicted  in 
Figure  2-3b. 

To  complement  the  quantitative  stay  time  data  employed  in  con¬ 
struction  of  the  Attrition  Rate  Model,  qualitative  behavior  profiles  and 
motivations  for  shelter  exit  as  a  function  of  time  were  also  derived  from 
the  literature  surveyed.  These  are  summarized  in  Table  2-2  for  the  first 
seven  days  following  the  attack  Generally,  post  attack  behavior  falls  in 
the  range  considered  normal.  As  shown  in  the  table,  reasons  for  shelter 
exit  include:  physical  and  physiological  effects,  such  as  shelter  damage 
and  lack  of  food,  water,  plumbing,  etc,;  lack  of  communication  and  direction; 
concern  for  relatives  and  friends;  boredom;  the  desire  to  explore  the  ex¬ 
ternal  environment;  the  perception  of  the  danger  being  over;  and  the  need 
for  action. 

2.4  DISCUSSION 

The  Attrition  Rate  Model,  as  denoted  in  this  study,  is  the  re¬ 
sponse  depicted  in  Figure  2-2.  Based  on  the  systematic  review  of  behavioral 
science  source  literature  and  the  development  of  a  database  of  quantitative 
stay  time  estimates,  the  model  indicates  that  the  expected  behavioral  re¬ 
sponse  of  a  sheltered  population  is  one  of  constant  “attrition"  from  the 
shelter  environment.  This  is  to  be  contrasted  with  traditional  assumptions 
in  which  the  population  is  taken  to  be  either  completely  in  or  completely 
out  of  their  protective  shelters  at  any  one  time. 

Because  of  its  foundation  in  a  quantitative,  empirical  database, 
the  Attrition  Rate  Model  lends  more  justification  to  its  use  as  a  fallout 
shelter  stay  time  assumption  then  do  currently  employed  assumptions  which 
vary  widely  in  their  definition,  and  which  represent  merely  reasonable  or 
convenient  model  imputs.  The  model  provides  a  significant  refinement  of 
fallout  shelter  stay  time  representation,  and  contributes  to  the  reduction 
in  variability  of  the  total  number  of  simulation  fallout  injuries  and 
fatalities  characteristic  of  the  range  of  stay  time  assumption  currently 
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Table  2-2.  Behavioral  profiles  and  reasons  for  shelter  exit. 
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- 

- 
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employed.  Although  the  model  is  based  on  criteria  which  attempt  to 
approximate  as  closely  as  possible  a  post-attack  environment,  and  upon  a 
database  which  was  specifically  sought  to  reflect  behavioral  response 
due  to  real  external  threats  and  a  spectrum  of  shelter  types,  applica¬ 
tion  of  the  model  to  situations  exceeding  the  characteristics  of  its 
database  must  be  made  with  care. 

For  example,  the  Attrition  Rate  Model  is  based  on  data  almost 
entirely  representative  of  American  response  to  disaster.  Thus,  appli¬ 
cation  of  the  model  to,  say,  the  Soviet  population  entails  a  cross-cultural 
extropolation  and  assumption.  One  may  estimate  that  the  commonality  of 
human  response  to  disaster  is  such  that  all  cultural  response  will  lie 
somewhere  within  the  Best  and  Worst  cases  displayed  by  the  model.  However, 
in  the  absence  of  additional  data,  this  must  remain  an  assumption.  In 
addition,  the  model  displays  a  clear  distinction  between  the  response  of 
well-trained,  well-equipped,  and  well-managed  shelterees  and  the  response 
of  those  ill-suited  in  these  areas.  Thus,  any  determined  or  overt 
program  which  tends  to  significantly  accentuate  these  positive  attributes 
prior  to  sheltering  may  condition  behavior  atypical  of  that  displayed  by 
the  database  employed.  However,  the  methodology  employed  in  constructing 
the  model  is  largely  independent  of  the  database  compiled;  just  the 
question  of  data  availability  remains.  In  constructing  the  Attrition 
Rate  Model ,  only  American  disaster  data  were  available. 

Because  the  model  indicates  a  continuous  attrition  of  the  shel¬ 
terees  from  the  shelter  environment,  application  of  model  results  within 
strategic  simulations  will  require  some  process  of  discretization.  In  ad¬ 
dition,  among  the  eight  categories  of  data  identified,  only  responses 
characteristic  of  bi-polar  parameters  were  developed.  For  example,  for  the 
case  of  training,  responses  reflecting  either  "presence  of  training"  or 
"absence  of  training"  are  represented.  Thus,  the  model  does  not  charac¬ 
terize  response  as  a  function  of  degrees  of  training,  or  levels  of  communi¬ 
cation,  or  completeness  of  warning,  etc.  However,  aside  from  issues  of 
model  domain  application  of  the  model  itself  within  currently  employed 
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strategic  simulations  is  envisioned  as  a  straightforward  exercise, 
involving  at  most  computational  and  integration  effort. 

Consequently,  the  Attrition  Rate  Model  is  presented  in  this  study 
as  an  "experimental  tool"  for  use  by  strategic  analysts  in  evaluating 
fallout  casualties  in  nuclear  exchange  simulations.  For  this  purpose,  the 
remainder  of  the  report  focuses  on  the  detailed  use  of  the  model  in  simu¬ 
lations  exhibiting  various  levels  of  detail.  In  particular,  two  appli¬ 
cations  are  presented:  model  use  in  a  "pencil-and-paper"  study  of  fallout 
casualties  in  Tbilisi,  Soviet  Georgia,  and  model  use  in  SIDAC. 
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SECTION  3 

APPLICATIONS  OF  THE  ATTRITION  RATE  MODEL 


The  usefulness  of  the  Attrition  Rate  Model  of  fallout  shelter 
stay  times  developed  in  the  previous  section  as  a  tool  for  the  investiga¬ 
tion  of  fallout  casualties  in  strategic  simulations  is  illustrated  in 
this  section.  First,  a  generic  guide  for  the  application  of  the  model 
is  presented.  This  guide  provides  a  framework  within  which  strategic 
studies  at  various  levels  of  detail  can  make  use  of  the  model.  Appli¬ 
cation  of  the  model,  and  use  of  the  guide,  is  then  presented  by  means  of 
two  examples:  first,  a  1 penci 1-and-paper'  study  of  fallout  casualties 
in  Tbilisi,  Soviet  Georgia,  and  second,  a  study  of  casualties  using  the 
SIDAC  simulation  model. 

3.1  A  GENERIC  GUIDE  TO  THE  USE  OF  THE  ATTRITION  RATE  MODEL 

This  section  describes,  in  general  terms,  procedures  for  the 
implementation  of  the  Attrition  Rate  Model  in  the  calculation  of  strate¬ 
gic  fallout  casualties.  The  steps  as  described  below  were  performed  in 
the  Tbilisi  study.  In  the  SIDAC  study,  only  the  last  step  was  specif¬ 
ically  performed;  the  remainder  are  performed  automatically  in  the  SIDAC 
simulation  itself. 

The  necessary  steps  are  as  follows: 

(1)  Determine  the  expected  distribution  of  the  post-attack 

population  among  fallout  shelters. 

(2)  Generate  the  weapon  laydown  and  develop  the  resulting 

blast  and  fallout  contours. 

(3)  Identify  those  shelters  destroyed  by  blast,  and  the 

distribution  of  surviving  shelters  among  fallout  radiation 
field  intensities. 

(4)  Employing  the  Attrition  Rate  Model,  determine  the  radia¬ 

tion  dose  received  by  each  individual  as  a  function  of 
his  shelter  protection  factor  and  the  intensity  of  the 
local  fal lout  field. 

(5)  Compute  casualties  based  on  dose  received. 
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These  guidelines  are  purposefully  general  in  order  that  they 
may  be  applicable  over  a  wide  range  of  simulation  detail.  In  the  examples 
to  follow,  use  of  the  model  is  keyed  to  investigating  the  sensitivity  of 
the  total  number  of  fallout  casualties  to  variations  in  stay  time  assump¬ 
tion. 

3.2  FIRST  MODEL  APPLICATION:  TBILISI,  SOVIET  GEORGIA 

The  study  of  a  simulated  laydown  pattern  on  Tbilisi,  Soviet 
Georgia,  was  a  preliminary  pencil-and-paper  study  to  demonstrate  the  use 
of  the  model.  The  goal  of  the  study  was  to  compare  results  for  two 
traditional  stay  time  assumptions  with  the  results  for  the  new  Attrition 
Rate  Model.  The  two  traditional  assumptions  investigated  were  the  Seven 
Day  assumption  (7  days  in,  14  days  2/3  in)  and  the  Two  Day  assumption  (2 
days  in,  3  days  2/3  in).  These  were  compared  with  the  Best  and  Worst 
Cases  from  the  Attrition  Rate  Model.  The  Best  and  Worst  Case  attrition 
rates  are  shown  in  Figure  3-1.  Only  a  summary  of  the  study  is  given 
here;  further  details  may  be  found  in  (3). 

3.2.1  Assumptions 

In  order  to  perform  the  comparison,  all  assumptions  of  the 

model  except  that  of  shelter  stay  times  were  held  constant.  For  the 

hypothetical  attack,  the  laydown  consisted  of  twelve  strategic  RVs  tar¬ 
geted  on  significant  industrial  or  military  targets  in  and  around  the 

city.  Six  RVs  were  employed  as  ground  bursts;  the  remaining  six  RVs 

were  air  bursts.  A  simultaneous  burst  laydown  was  assumed. 

The  fallout  was  caused  by  a  20  knot  wind  blowing  northwest,  a 
direction  characteristic  of  the  area  between  November  and  April.  Fallout 
was  assumed  to  arrive  at  1.5  hours.  (This  is  reasonable,  as  we  only 
considered  the  casualties  within  the  city;  downwind  rural  casualties 
were  not  investigated). 

The  distribution  of  population  among  shelters  of  various  PFs 
and  the  shelter  blast  vulnerabilities  were  taken  from  an  SRI  study  (4). 
The  shelters  were  assumed  to  be  uniformly  distributed  within  the  city 
boundaries.  The  outside  protection  factor  was  assumed  to  be  PF=2. 
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Days  Since  Sheltering 

Best  Case 
y  =  2.5  +  3.5t 

Trained 

Warned 

Food  and  Water  Generally 
Available 
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Figure  3-1.  The  best  and  worst  case  attrition  rate  models. 
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To  estimate  shelter  destruction  by  blast  damage,  the  DIA  blast 
vulnerability  methodology  (Physical  Vulnerabilties  Handbook  (5))  was  used. 
Fallout  radiation  fields  were  modeled  using  the  EM-1  fallout  methodology  (6). 
Fallout  was  assumed  to  decay  as  t  '  ,  with  t  in  hours. 

The  number  of  casualties  was  equal  to- 

N=  (1/4)' (individuals  receiving  between  250R  and  450R) 

+(3/4)‘ (those  receiving  between  450R  and  650R) 

+(1.0)' (those  receiving  over  650R) 

All  doses  were  received  within  ten  days,  and  in  most  cases  the  majority  of 
the  dose  was  received  within  four  days,  so  no  biological  repair  factor  was 
needed  or  used. 

For  ease  of  computation,  attrition  from  fallout  shelters  was 
assumed  to  take  place  at  24  hour  intervals,  rather  than  continuously.  This 
introduces  no  significant  change  in  the  results. 

3.2.2  Approach 

To  estimate  the  variation  in  sheltered  population  fallout 
casualties  in  Tbilisi  for  each  of  the  stay  time  assumptions,  a  two  part 
approach  was  employed. 

In  the  first  part,  the  number  of  casualties  per  day  was  computed 
for  shelters  with  PFs  of  10,  50,  100,  200  and  500  located  in  fallout  fields 
ranging  in  intensity  (at  one  hour)  from  300  to  10,000  rads  per  hour. 
Casualties  were  computed  for  the  four  different  stay  time  assumptions:  Seven 
Day,  Two  Day  and  the  Best  and  Worst  Cases  from  the  Attrition  Rate  Model. 

In  the  second  part,  a  hypothetical  attack  on  Tbilisi  was  consider¬ 
ed.  A  laydown  pattern  was  generated  as  discussed  above.  The  resulting  blast 
damage  to  shelter  was  computed  using  the  Physical  Vulnerabilities  Handbook. 
Figure  3-2  shows  the  city  boundaries  and  regions  of  shelter  destruction. 
Fallout  contours  were  constructed  using  EM-1  methodology.  In  Figure  3-3  the 
radiation  field  resulting  from  the  laydown  is  presented.  The  fraction  of 
surviving  (undamaged)  shelters  in  each  radiation  field  was  computed  by 
overlaying  each  part  of  Figure  3-2  with  Figure  3-3  and  computing  the  area 
contained  in  each  field.  The  resulting  number  of  daily  fallout  casualties 
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was  then  determined  by  estimating  the  percent  of  undestroyed  shelters  in 
fallout  fields  of  intensity  300  to  10,000  R/hr,  and  employing  the  results  of 
part  one. 

3.2.3  Results 

Part  one  resulted  in  four  graphs  for  each  protection  factor:  one 
graph  for  each  of  the  four  possible  stay  time  models.  As  an  example,  Figure 
3-4  presents  the  results  for  PF=200.  (Complete  results  may  be  found  in  (3)). 

Note  that  the  two  versions  of  the  Attrition  Rate  Model  yield 
relatively  smooth  curves.  This  is  because  of  the  continuous  outflow  of 
shelterees  under  this  model.  The  traditional  assumptions,  on  the  other  hand, 
are  much  more  jagged.  In  fact,  by  definition  these  curves  can  only  take 
values  of  0,  25,  75  and  100  percent. 

The  results  of  part  two  are  summarized  in  a  series  of  tables. 
Table  3-1  provides  a  summary  of  the  distribution  of  the  sheltered  population 
both  before  and  after  the  blast.  In  Table  3-2  the  population  which  survived 
the  blast  is  distributed  by  percentage  among  the  various  protection  factors 
and  radiation  fields.  Figure  3-5  and  Table  3-3  depict  the  final  casualty 
figures.  (Note;  In  Figure  3-5  and  Table  3-3,  a  casualty  is  defined  as  one  who 
has  received  a  fatal  dose,  not  necessarily  one  who  has  already  died.  The 
actual  times  of  death  might  be  later  than  depicted  in  the  graphics). 

3.2.4  Discussion 

A  number  of  conclusions  may  be  drawn  from  Figure  3-4  and  similar 
diagrams  for  other  protection  factors.  Using  the  Two  Day  and  Seven  Day 
assumptions,  no  casualties  result  from  a  300  R/hr  or  less  radiation  field. 
Under  the  Attrition  Rate  assumption,  those  people  who  leave  shelter  within 
the  first  few  hours  were  casualties.  (These  individuals  were  assumed  to  be 
sheltered  from  the  blast,  but  to  leave  shelter  before  the  fallout  arrived.) 
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Figure  3-4.  Daily  casualties  for  the  first  ten  post-attack  days  as  a  function  of 

stay  assumption  and  fallout  radiation  intensity  (in  rads/hr  at  one  hour) 
for  shelters  with  PF  =  200. 
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Figure  3-5.  Daily  percentage  of  sheltered  Tbilisi  population  becoming 
casualties  within  the  first  ten  days  for  each  of  the  three 
shelter  stay  time  assumptions. 


3-3.  Fallout  casualties  in  Tbilisi  (percent). 
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The  major  difference  between  the  casualty  curves  for  the  tradi- 
nal  and  Attrition  Rate  assumptions  occurs  for  high  PF  shelters.  The  tradi¬ 
tional  assumptions  result  in  100%  casualties  for  those  in  high  radiation 
fields,  and  100%  survival  for  those  in  low  radiation  fields.  The  Attrition 
Rate  assumptions  more  realistically  allow  other  than  "all-or-nothing" 
choices.  Casualties  are  higher  in  high  fields  than  in  low  ones,  but  those 
people  who  remain  in  high  PF  shelters  have  a  chance  to  survive  in  even  the 
worst  radiation  fields. 

As  indicated  in  Figure  3-5,  the  percentage  of  post-attack  casual¬ 
ties  associated  with  the  Attrition  Rate  stay  time  assumption  is  between  1.5 
to  2  times  greater  than  the  percentage  expected  for.  the  Seven  Day  assump¬ 
tion,  and  between  1/8  to  1/2  times  less  than  casualties  expected  based  on 
the  Two  Day  assumption  (after  the  fourth  day).  However,  the  Attrition  Rate 
assumption  is  based  on  a  survey  of  human  behavioral  response  across  a 
spectrum  of  shelter  types,  whereas  the  Two  Day  and  Seven  Day  cases  assume, 
respectively,  shelters  characterized  by  inadequate  and  generally  adequate 
water  availability  and  plumbing  facilites.  Thus,  any  comparison  of  these 
stay  time. assumptions  must  include  the  expected  distribution  of  Two  Daylike 
and  Seven  Daylike  shelters  within  the  country  under  attack. 

According  to  the  SRI  study  "Fallout  Protection  Factors  for  Various 
Countries"  (4),  shelters  in  the  U.S.  and  Soviet  Union  are  overwhelmingly 
of  the  Seven  Day  type.  Two  Daylike  shelters  are  more  prevalent  in  lesser 
developed  countries  such  as  Spain,  Turkey,  Romania,  Hungary,  etc.  As  a 
result,  this  analysis  indicates  that  the  "standard"  stay  time  assumption 
for  the  U.S.  and  Soviet  Union  (i.e.,  Seven  Day)  probably  underestimates 
post-attack  casualties  by  between  five  and  fifteen  percentage  points.  In 
addition,  as  shown  by  the  Worst  Case  and  Best  Case  results  of  Figure  3-5, 
adequate  training,  warning,  and  general  preparedness  would  reduce  expected 
casualties  between  five  to  ten  percentage  points  in  Tbilisi.  For  a  city 
of  this  size,  this  would  represent  approximately  60,000  civilians. 
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In  some  cases,  an  analyst  would  use  both  the  Two  Day  and  Seven 
Day  stay  times  to  investigate  the  variability  of  casualties  as  a  function 
of  stay  time  assumptions.  There  is  a  250%  difference  between  the  Two  Day 
and  Seven  Day  results  after  ten  days.  The  Best  and  Worst  Cases,  however, 
seem  to  reduce  this  variability  to  approximately  50%. 

It  should  be  noted  that  the  calculations  performed  during  this 
study  have  shown  that  the  daily  percentage  of  casualties  expected  for  the 
Attrition  Rate  and  Seven  Day  assumptions  is  a  function  of  the  shelter 
protection  factor  and  the  intensity  of  the  local  radiation  field.  That 
is,  the  conclusion  cannot  be  made  that  for  all  PFs  and  fallout  inten¬ 
sities,  for  example,  the  Attrition  assumption  will  always  generate  more 
casualties  than  the  Seven  Day  assumption.  Consequently,  the  variation 
in  casualties  for  each  assumption  is  expected  to  be  scenario  dependent, 
requiring  the  inclusion  of  estimated  shelter  distributions  within  con¬ 
structed  fallout  fields,  and  consideration  of  population  allocations  to 
those  shelters. 

In  summary,  the  exemplative  analysis  of  an  attack  on  Tbilisi 
conducted  for  this  study  has  shown: 

(1)  The  number  of  post-attack  casualties  is  expected  to  be 
scenario  dependent; 

(2)  The  inclusion  of  realistic  behavioral  assumptions  in  stay 
times  (Attrition  Rate  assumption)  generated  1.5  to  2 
times  more  casualties  in  Tbilisi  than  the  "traditional" 
assumption  (Seven  Day),  and; 

(3)  Adequate  shelter  preparedness  generated  between  1/4  and 
1/3  fewer  casualties  than  inadequate  shelter  preparedness 
in  Tbilisi. 

(4)  The  Attrition  Rate  Model  reduces  the  range  of  variability 
associated  with  stay  time  assumptions. 
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3.3  SECOND  MODEL  APPLICATION:  SIDAC 

3.3.1  The  SIDAC  Model 

SIDAC  is  a  computerized  analytical  model  designed  to  provide 

nuclear  damage  anlaysis  information  for  both  Red  and  Blue  resource 

monitoring.  It  is  a  one-sided  model  that  simulates  land,  air,  and  sea 

forces,  as  well  as  civilians  and  paramilitary.  It  can  consider  weapons 
or  weapons  systems  individually  and  the  modularity  of  its  design  allows 
the  user  to  aggregate  up  to  any  level  he  wishes,  depending  upon  his 
specific  requirements.  The  model  uses  a  mixture  of  deterministic  and 

stochastic  elements.  Probability  is  used  as  the  primary  solution  tech¬ 
nique  for  prompt  damage  by  means  of  the  methodology  developed  by  the 

Physical  Vulnerability  (PV)  Division  of  the  United  States  Air  Force 

Intelligence.  Delayed  radiation  effects  are  estimated  by  means  of  the 
methodology  developed  by  the  Weapons  Systems  Evaluation  Group  (WSEG). 

SIDAC  was  developed  by  the  Command  and  Control  Technical 
Center  (CCTC)  of  the  Defense  Communications  Agency.  It  is  used  by  the 
Studies,  Analysis,  and  Gaming  Agency  (SAGA)  under  the  aegis  of  the  Joint 
Chiefs  of  Staff. 

The  structure  of  SIDAC  is  summarized  in  Figure  3-6.  Inputs  to 
SIDAC  consist  of  the  strike  tape  (containing  the  weapon  laydown  and 
related  weapon  information),  wind  and  weather  conditions,  and  the  data 

base,  containing  target  and  population  information. 

SIDAC  uses  this  information  to  produce  an  Answer  File,  often 

designated  by  File  Code  25.  For  our  purposes,  we  may  consider  the 
Answer  File  to  have  one  "record"  (unit  of  output)  for  each  possible 
environment.  An  environment  may  be  considered  to  be  a  group  of  fallout 
shelters  (of  varying  PFs)  in  a  given  radiation  field  with  a  given  prob¬ 
ability  of  blast  damage. 

A  record  contains  basic  identification  data,  such  as  the 
country  and  geopolitical  region  of  the  information  and  whether  it  is  in 
an  urban  or  rural  environment.  Also  included  are  the  capacity  of  the 
shelter  (CAP),  the  probability  of  receiving  less  than  moderate  or  severe 
blast  damage  (MPROB  and  SPROB,  respectively),  and  the  maximum  cumulative 
biological  dose  (MAXDOS)  which  would  be  obtained  by  an  individual  in 
that  environment  with  no  fallout  shielding  (PF=1). 
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Figure  3-6.  SIDAC  structure  and  information  flow. 


3.3.2 


The  SIDAC  Post-Processor 

The  SIDAC  post-processor  performs  the  actual  computation  of  prompt 
and  fallout  casualties.  The  post-processor  as  originally  developed  by  CCTC 
consists  of  approximately  350  lines  of  FORTRAN  code.  A  listing  is  contained 
in  Appendix  C.  This  program  uses  the  SIDAC  Answer  File  as  input  and  produces 
casualty  estimates  for  the  scenario  under  consideration. 

As  originally  configured,  the  post-processor  allowed  the  user  to 
input  a  shelter  stay  time  assumption  to  apply  to  the  entire  population.  This 
stay  time  is  described  in  two  parts:  the  period  from  time  0  to  time  ITT 
represents  the  period  of  100%  shelter  occupancy.  The  period  from  ITT  to  JTT 
represents  the  period  of  partial  shelter  occupancy.  The  fraction  of  time  the 
shelter  is  occupied  is  designated  XMULT.  After  time  JTT,  radiation  exposure 
was  assumed  to  drop  to  zero.  This  is  due  either  to  evacuation  of  the  area  at 
risk  or  actual  decay  of  radiation  levels  to  below  noticable  limits. 

An  example  will  serve  to  clarify  this.  Assume  we  want  to 
investigate  the  effects  of  a  7  day  shelter  stay  time  followed  by  14  days  of 
two-thirds  sheltering.  After  21  days,  evacuation  takes  place.  All  times 
must  be  in  hours.  Therefore,  we  set  ITT  =  7  x  24  =  168  hours,  JTT  =  21  x  24  = 
504  hours  and  XMULT  =  .667.  With  these  inputs,  the  program  would  calculate 
the  resulting  casualties.  Further  details  are  contained  in  Appendix  B. 

3.3.3  The  Modified  Post-Processor 

The  post-processor  as  described  above  was  modified  to  allow  a 
linear  attrition  rate  of  y  =  at  +  b,  as  described  in  the  Attrition  Rate  Model, 
with  t  =  time  in  hours  and  y  =  proportion  leaving  shelter.  The  modifications 
will  be  described  in  brief  here  and  more  fully  in  Appendix  D. 

The  Attrition  Rate  Model  calls  for  a  continuous  flow  of  shelterees 
from  the  shelter.  For  computational  purposes,  this  was  approximated  by  a 
discrete  model  with  one  exodus  from  the  shelters  every  24  hours.  For 
example,  consider  the  linear  attrition  equation  y  =  . 00 i 5 t  +  .025.  In  this 
case,  2.5%  of  the  shelterees  leave  immediately,  and  an  additional  .0015  x  24 
=  3.6%  of  the  shelterees  leave  every  24  hours. 
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To  implement  this  model,  a  loop  was  inserted  in  the  original 
program.  ITT  was  started  at  0  and  incremented  by  24  hours  for  each  new  day's 
population  leaving  shelter.  JTT  was  held  fixed  at  504  hours  (21  days)  as  it 
always  was  in  the  original  version.  The  Attrition  Rate  Model  makes  no 
allowance  for  partial  sheltering,  so  XMULT  =  0. 

For  each  value  of  ITT,  casualty  computations  were  made  as  in  the 
original  post-processor.  In  this  case,  however,  the  casualty  figures  were 
multiplied  by  that  fractions  of  the  population  actually  leaving  shelter  at 
time  ITT.  The  casualties  for  groups  leaving  shelter  on  each  day  were  summed 
to  give  the  total  casualty  figures. 

Again,  an  example  will  serve  to  clarify  matters.  Consider,  for 

25 

simplicity  the  attrition  rate  equation  y  =  t  +  .10,  with  t  in  hours. 
This  says  that  10%  of  the  population  leaves  shelter  immediately,  and  an 
additional  25%  leaves  every  24  hours  until  the  shelters  are  empty.  After  72 
hours,  85%  have  left.  On  day  4  (96  hours)  the  remaining  15%  leave.  We 
further  assume  that  the  radiation  field  is  such  that  100%  of  those  leaving 
immediately  die,  as  do  60%  of  those  leaving  after  1  day,  40%  of  those  leaving 
after  2  days,  20%  of  those  leaving  after  3  days,  and  10%  of  those  leaving 
after  4  days.  Consider  a  sheltered  population  of  1000  people.  Table  3-4 
summarizes  the  calculations.  Out  of  1000  people,  415  became  fatalities. 

The  modified  post-processor  calculates  casualties  for  an  arbitrary 
attrition  rate  y  =  at  +  b.  The  two  required  inputs  are  the  parameters  a  and 
b.  The  modifications  to  the  original  code  consist  of  approximately  25  lines 
out  of  350.  The  modified  post-processor  appears  in  Appendix  0. 

3.3.4  Method  of  Determining  Casualties 

The  original  and  modified  post-processors  use  almost  identical 
methods  to  compute  casualties.  In  fact,  the  only  difference  is  that  the 
modified  version  multiplies  casualties  from  a  given  day's  attrition  by  the 
proportion  leaving  on  that  day  and  then  sums  across  days;  the  original 
assumes  all  attrition  occurs  at  once,  and  so  only  has  one  group  to  consider. 
Therefore,  we  only  consider  the  method  in  the  original  post-processor. 
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Table  3-4. 

Casualty 
y  =  0.25t 

calculation  for  hypothetical  attrition  rate 
+  0.10 

DAY 

PEOPLE 

LEAVING 

FATALITY 
x  FRACTION 

NUMBER  OF 

FATALITIES 

0 

100 

1.00 

100 

1 

250 

.60 

150 

2 

250 

.40 

100 

3 

250 

.20 

50 

4 

150 

.10 

15 

TOTAL 

1000 

415 

I 


46 


All  injury  and  fatality  calculations  are  made  for  three  time 
periods:  7  days,  30  days,  and  180  days  after  the  blast.  Unlike  the 
Tbilisi  study  discussed  above,  this  program  counted  a  fatality  only  when 
it  actually  occurred,  not  when  the  lethal  dose  was  received. 

Recall  that  each  record  in  the  Answer  File  describes  one 
outside  radiation  field  and  level  of  blast  damage.  The  population  des¬ 
cribed  in  this  record  is  assumed  to  be  distributed  among  shelters  of 
various  PFs,  as  shown  in  Table  3-5. 

The  overall  structure  of  casualty  determination  is: 

(1)  Compute  the  proportion  of  prompt  fatalities  and  injuries 
based  on  probabilities  of  damage. 

(2)  For  each  radiation  field  and  protection  factor,  compute 
the  proportion  of  fallout  fatalities  and  injuries. 

(3)  Multiply  number  of  shelterees  by  proportion  of  prompt 

fatal ities  to  obtain  the  number  of  prompt  fatal i ties . 

(4)  Multiply  number  of  remaining  shelterees  (not  killed  by 

prompt  effects)  by  proportion  of  fallout  casualties  to 
obtain  the  number  of  fa! 1  out  fatal ities. 

(5)  Multiply  number  of  still  remai ni ng  shelterees  by  propor¬ 
tion  of  prompt  injuries  to  obtain  the  number  of  prompt 

injuries. 

(6)  After  subtracting  the  number  of  fatalities  and  prompt 
injuries,  multiply  the  number  of  remaining  snelterees  by 
the  proportion  of  fallout  injuries  to  obtain  the  number 
of  fal lout  i njuries .  Note:  Steps  5  and  6  imply  that 
joint  prompt  and  fallout  injuries  are  counted  simply  as 
prompt  injuries. 

(7)  The  uninjured,  healthy  population  is  the  remaining  popu¬ 
lation  after  subtracting  all  casualties. 

(8)  This  procedure  is  performed  for  each  data  record,  and  the 
results  are  summed. 

This  procedure  is  described  in  more  detail  in  Appendix  B. 
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Table  3-5.  Distribution  of  shelterees  for  urban  and  rural  environments. 


URBAN  RURAL 


PF 

PERCENT 

PF 

PERCENT 

800 

14 

40 

2 

250 

7 

35 

5 

150 

7 

30 

15 

100 

14 

25 

8 

40 

28 

20 

22 

20 

17 

15 

6 

10 

13 

10 
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*  These  shelter  distributions  may  be  explicitly  overruled  by  input 
data,  but  this  was  never  done  during  our  sample  runs. 


3.3.5 


The  SIDAC  Run 

CCTC  provided  a  SIDAC  Answer  File  for  analysis.  The  scenario  was 
based  on  a  Soviet  attack  on  the  U.  S.  For  classification  reasons,  details 
of  the  SIDAC  scenario  are  omitted.  Because  of  this,  the  actual  casualty 
figures  should  not  be  considered  as  representative  of  all  SIDAC  runs. 
However,  the  relative  spread  of  results  is  significant. 

Fallout  shelter  stay  times  were  varied  to  investigate  the 
sensitivity  of  casualty  figures  to  stay  times.  Nine  different  assumptions 
were  used: 

(1)  3/4  day  fully  sheltered,  20k  days  2/3  sheltered 

(2)  2  days  fully  sheltered,  19  days  2/3  sheltered 

(3)  3  days  fully  sheltered,  18  days  2/3  sheltered 

(4)  5  days  fully  sheltered,  16  days  2/3  sheltered 

(5)  7  days  fully  sheltered,  15  days  2/3  sheltered 

(6)  14  days  fully  sheltered,  7  days  2/3  sheltered 

(7)  21  days  fully  sheltered 

(8)  Attrition  rate  Best  Case,  y  =  .0014t  +  .025,  no  partial 

sheltering  J/ 

(9)  Attrition  rate  Worst  Case,  y  =  .0025t  +  .055,  no  partial 

sheltering. 1/ 

For  all  nine  cases,  radiation  exposure  is  assumed  to  end  at  21  days  due 
to  evacuation  or  the  decay  of  radiation  intensity  to  insignificant  levels. 
Note  that  this  means  No.  7  is  equivalent  to  indefinite  sheltering. 

It  is  instructive  to  consider  alternative  stay  times  that 
provide  identical  effective  protection  factors  for  the  21-day  period. 
Some  of  these  are  shown  in  Table  3-6. 

3.3.6  Results 

The  SIDAC  data  base  assumes  a  total  U.S.  population  of  214.6 
million.  Of  these,  131.4  million  are  urban,  while  83.2  million  are  rural. 
The  post-processor  provides  separate  casualty  figures  for  the  urban  and 
rural  populations. 


1.  For  Nos.  8  and  9,  t  is  in  hours. 
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Table  3-6.  Effective  protection  factors  and  equivalent  exit  days. 


PF 

JO 

100 

250 

800 

SIDAC  STAY  TIME 
ASSUMPTION 

PFe 

EXIT 

DAY 

PFe 

EXIT 

DAY 

PFe 

EXIT 

DAY 

PFe 

EXIT 

DAY 

2  days  fully 
sheltered,  19 
2/3  sheltered 

days 

6.6 

8.2 

15 

8.6 

16 

8.5 

17 

8.6 

3  days  fully 
sheltered,  18 
2/3  sheltered 

days 

7.2 

10.2 

18 

10.2 

20 

10.2 

21 

10.1 

5  days  fully 
sheltered,  16 
2/3  sheltered 

days 

7.8 

12.3 

23 

12.3 

27 

12.4 

29 

12.3 

7  days  fully 
sheltered,  14 
2/3  sheltered 

days 

8.3 

14.2 

29 

14.1 

35 

14.1 

39 

14.2 

DEFINITIONS: 

PFg  =  Effective  protection  factor 

Equivalent  exit  day  =  the  day  such  that  complete  sheltering  up  to 
that  day  and  no  sheltering  after  it  gives  the  same  PF  as  the 
corresponding  SIDAC  stay  time  assumption. 

ASSUMPTIONS: 

Fallout  arrives  at  1^  hours 
Outside  PF  *  1/.65  =  1.538 
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The  actual  post-processor  output  is  provided  in  Appendix  E. 
This  section  contains  a  summary  and  analysis  of  the  results. 

Prompt  fatality  calculations  are  independent  of  subsequent 
fall  out  shelter  stay  times.  Therfore,  every  set  of  results  has  the 

same  number  of  prompt  fatalities.  For  this  scenario,  92.1  million  of 
the  urban  population  (70.1%)  and  7.5  million  of  the  rural  population 
(9.0%)  were  prompt  fatalities. 

Prompt  injuries  are  slightly  dependent  on  fallout  sheltering. 

(A  prompt  injury  can  become  a  fallout  fatality,  and  this  can  happen  at 
different  times  for  different  levels  of  sheltering.)  However,  the 
figures  are  relatively  constant  over  time  and  for  each  scenario.  For 
comparative  purposes,  figures  of  23.0  million  urban  prompt  injuries 
(17.5%)  and  11.5  million  rural  prompt  casualties  (13.8%)  were  used. 

Actual  results  differed  from  these  by  no  more  than  two  or  three  percent¬ 
age  points.  The  actual  results  are  available  in  Appendix  E. 

The  post-processor  provides  casualty  figures  for  three  differ¬ 
ent  points  in  time:  7  days,  30  days,  and  180  days  after  the  attack. 
Fallout  fatalities  and  injuries  are  assessed  on  the  basis  of  maximum 
biological  dose  (MBD)  received.  In  cases  of  shorter  shelter  stay  times 
(less  than  5  days)  this  M8D  is  received  before  the  seventh  day,  and  so 
the  7  day  casualty  figures  are  accurate.  However,  for  longer  stay 
times,  this  MBD  is  not  received  until  sometime  after  seven  days  have 
elapsed,  so  the  7  day  fallout  casualty  figures  are  inaccurate. 

Because  of  the  high  percentage  of  prompt  casualties,  only  a 

small  fraction  of  the  total  population  was  at  risk  (i.e.,  alive  after 
the  blast)  for  fallout  casualty  calculations.  In  addition  this  fraction 
was  different  for  urban  and  rural  populations.  Therefore,  in  presenting 
fallout  fatalities,  the  percentages  were  normalized  to  the  population  at 
risk  by  dividing  by  the  fraction  of  population  which  were  not  prompt 
fatalities.  Similarly,  for  fallout  casualties,  the  percentages  were 
normalized  by  dividing  by  the  fraction  not  prompt  fatalities  or  injuries. 


With  these  points  in  mind,  consider  the  casualty  figures 
presented  in  Figure  3-7  and  3-8.  Figure  3-7  presents  fallout  fatalities 
for  the  nine  stay  times  and  for  both  the  urban  and  rural  populations 
(left  and  right  graphs,  respectively).  Figure  3-8  presents  similar 
information  for  fallout  injuries. 

Notice  that  the  error  in  the  7  day  figures  discussed  above  is 
readily  apparent.  If  there  were  no  error,  the  7  day,  14  day,  and  21  day 
stay  time  figures  would  be  identical.  The  differences  demonstrate  that 
the  MBD  is  not  always  received  by  the  seventh  day. 

These  graphs  effectively  illustrate  the  major  goal  in  developing 
the  Attrition  Rate  Model:  reduction  in  variability  associated  with 
sheltering  estimates.  At  180  days,  there  is  a  500%  difference  in  urban 
fatalities  and  a  250%  difference  in  rural  casualties  between  the  21  day 
full  sheltering  assumption  and  the  3/4  day  full,  1 9-s  day  2/3  sheltering 
assumption.  The  difference  between  the  Best  and  Worst  Case  is  roughly 
50%  for  urban  fatalities  and  25%  for  rural  fatalities.  Similar  relative 
differences,  although  much  smaller  in  actual  numbers,  exist  for  the 
injury  graph. 

The  Worst  and  Best  Cases  are  approximated  by  the  2  day  fully 
sheltered,  19  day  partial  sheltered  and  the  3  day  sheltered,  18  day 
partial  sheltered  assumptions.  However,  note  that  from  Table  3-6,  the  2 
day  sheltered,  19  day  partially  sheltered  stay  time  provides  the  same 
protection  as  an  8.6  day  fully  sheltered  stay,  followed  by  12.4  days  out 
in  the  local  radiation  field  (for  a  total  of  21  days)  before  evacuation. 
The  3  day  sheltered,  18  day  partially  sheltered  stay  time  is  equivalent 
to  10.2  days  in  a  shelter,  followed  by  10.8  days  in  the  local  radiation 
field  before  evacuation.  These  numbers  indicate  that  care  must  be 
exercised  in  attempting  to  compare  results  of  the  Attrition  Rate  Model 
with  one  "equivalent"  stay  time;  there  are  many  combinations  of  stay 
times  which  provide  equivalent  protection. 
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URBAN  POPULATION: 

131.4  MIUJON 

92.1  MILLION  (781NI  PROMPT  FATALITIES 
23  a  MILLION  limi  PROMPT  INJURIES 


RURAL  POPULATION: 

812  MILLION 

7.5  MILLION  ISPSI  PROMPT  FATALITIES 

11.5  MILLION  I11ISI  PROMPT  INJURIES 


14  DAY  SHELTERED. 

28 114  DAT  213  SHELTERED 


314  DAT  SHELTERED. 

28114  OAV  213  SHELTERED 


WORST  CASE 
2  DAYS  SHELTERED. 

IS  OATS  213  SHELTERED  ^1 

3  OATS  SHELTERED 
18  DAYS  2)3  SHELTEREO 
BEST  CASE 
5  OAYS  SHELTERED  ' 
18  OATS  213  SHELTEREO 

N 

7  OAYS  SHELTEREO 
14  OAYS  23  SHELTERED 

14  OAYS  SHELTEREO 
7  OAYS  213  SHELTEREO 
21  OAYS  SHELTEREO  * 


WORST  CASE 
2  DAYS  SHELTERED.  V 
19  OAYS  213  SHaTERB)  \ 

3  DAYS  SHELTEREO  v\ 
II  OAYS  213  SHELTERED  \  ' 
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Figure  3-8.  Fallout  injuries  for  urban  and  rural  populations. 
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There  is  a  significant  difference  in  the  computer  time  needed 
to  perform  these  runs.  The  runs  were  made  on  a  Honeywell  6000  series 

computer  running  under  the  GCOS  operating  system.  A  single  run  of  one 
of  the  Attrition  Rate  assumption  took  roughly  9  times  longer  than  a 
single  run  of  one  of  the  traditional  stay  time  assumptions.  However, 

recall  that  the  post-processor  was  designed  to  work  for  the  traditional 
assumptions,  and  "brute  force"  was  used  to  make  it  handle  the  Attrition 
Rate  assumptions.  A  program  properly  designed  specifically  for  the 
Attrition  Rate  assumptions  could  be  expected  to  improve  on  these  figures 
considerably. 

The  Best  and  Worst  Case  results  are  roughly  parallel  to  the 
more  traditional  stay  time  results  on  each  of  the  four  graphs.  They 
also  lie  roughly  in  the  same  range,  i.e.  ,  near  the  results  for  2  day 
and  3  day  stay  times.  However,  all  four  graphs  come  from  only  one  SIDAC 
scenario,  one  with  a  very  large  proportion  of  blast  casualties.  One  may 
hypothesize  that  these  relationships  hold  in  general  for  other  scenarios, 
or  that  the  3  day  fully  sheltered,  18  day  partially  sheltered  stay  time, 
for  example,  produces  the  same  number  of  casualties  as  the  Best  Case. 

To  date,  these  conjectures  remain  just  that,  and  further  work  is  needed 

to  substantiate  them. 

3.4  SUMMARY 

The  Attrition  Rate  Model  displays  three  advantages  as  a  tool 
for  the  investigation  of  fallout  casualties.  First,  it  has  a  basis  in  a 
real  data  base  analyzed  using  the  methods  of  the  behavioral  sciences. 
As  such,  it  provides  a  justification  for  the  use  of  model  stay  times;  a 
justification  absent  in  previous  stay  time  assumptions.  Second,  the 
Attrition  Rate  Model  reduces  the  variability  associated  with  a  variety 
of  shelter  stay  times.  Third,  the  methodology  used  to  derive  the  Attrition 
Rate  equations  can  be  applied  to  a  variety  of  behavioral  problems.  If 
appropriate  data  bases  are  available,  one  may  apply  these  methods  to 
develop  empirical  models  in  other  fields. 
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SECTION  4 

FURTHER  APPLICATIONS 


Among  the  initial  objectives  of  this  study  was  the  development  of  a 
tool  for  use  in  strategic  simulations  with  the  purpose  of  reducing  the 
variability  in  simulated  results  supported  by  the  "human  element".  The 
Attrition  Rate  Model  represents  such  a  tool  with  respect  to  fallout  shelter 
stay  times.  In  the  previous  section,  two  applications  of  the  model  were 
presented.  These  applications  focused  on  an  investigation  of  the  sensitivity 
of  fallout  casualties  to  stay  time  assumption.  In  this  chapter,  additional 
applications  of  the  model  are  suggested.  These  include  extensions  of  the 
type  of  research  exemplified  by  the  applications  in  Section  3,  as  well  as 
more  detailed  use  of  the  model. 

4. 1  IMPLEMENTATION  OF  THE  COMPLETE  SET  OF  ATTRITION  RATE  EQUATIONS 

The  first  application  is  to  modify  the  SIDAC  post-processor 
(Appendix  D)  to  accept  the  complete  set  of  24  different  Attrition  Rate 
equations.  These  equations  were  previously  presented  in  Table  2-4. 

There  are  two  distinct  parts  to  the  modification,  corresponding  to 
the  two  different  functional  forms  of  equations.  For  the  linear  equations  (y 
=  at  +  b)  almost  no  modification  is  needed.  The  program  is  already  designed 
to  handle  the  two  linear  equations  describing  the  Best  and  Worst  Cases.  The 
only  required  inputs  are  the  parameters  a  and  b  (designated  AA  and  BB  in  the 
modified  code,  Appendix  0).  By  inputting  the  proper  a  and  b,  the  current 
post-processor  will  handle  any  of  the  linear  equations. 

Of  the  24  equations,  5  take  the  simple  exponential  form  y  =  atb. 
This  is  slightly  more  difficult  to  handle  in  that  the  code  as  currently 
structured  cannot  handle  an  exponential  stay  time.  But  the  addition  of  this 
capability  is  fairly  straightforward. 

By  adding  this  capability,  variations  in  fallout  casualties  due  to 
the  range  of  response  in  each  of  the  eight  data  categories  summarized  in 
Table  2-4  could  be  investigated.  For  example,  the  impact  of  training  vs  no 
training,  or  good  shelter  management  vs  poor  shelter  management  on  population 
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survival  could  be  analyzed.  With  this  information,  the  analyst  could 
investigate  each  of  the  eight  variables,  rather  than  just  the  best  and  worst 
case  summaries. 

4.2  ADDITIONAL  SCENARIOS 

Tc  date,  only  one  SIDAC  scenario  has  been  run  as  described  in 
Section  3.3.  It  is  impossible  to  draw  general  conclusions  about  model 
behavior  on  the  basis  of  one  run.  Many  different  scenarios  must  be  consi¬ 
dered  and  the  variability  of  the  results  must  be  examined.  In  this  way  we  can 
investigate  various  hypothesis  about  shelter  stay  times  with  range  of 
scenarios  could  answer  questions  concerning: 

(1)  The  existence  of  simple  stay  time  assumptions  "equivalent"  to 
the  Best  and  Worst  Case  Attrition  Rate  Model; 

(2)  The  impact  of  the  Attrition  Rate  Model  assumption  for  various 
degrees  of  evacuation  posture; 

(3)  The  range  of  casualties  between  the  Best  and  Worst  case 
responses; 

(4)  The  variation  of  fallout  casualties  under  a  variety  of  weapon 
laydown  patterns;  and 

(5)  Other  hypothesis  concerning  the  model. 

It  is  a  simple  matter  to  implement  this  application.  Each  SIDAC  run  should 
be  made  in  the  usual  way,  generating  an  Answer  File.  This  Answer  File  is  then 
run  through  the  modified  post-processor  and  the  results  analyzed  as  in 
Section  3.3.  These  results  will  yield  general  rules  which  the  Attrition  Rate 
Model  follows. 

As  an  example,  consider  the  claim  that  a  3  days  sheltered,  18  days 
partially  sheltered  stay  time  yields  casualty  figures  approximately  equal  to 
the  Best  Case.  This  is  borne  out  in  the  one  scenario  that  has  been  studied. 
Should  this  prove  to  be  the  case  over  a  number  of  widely  varying  scenarios, 
we  would  consider  the  claim  validated.  If  this  claim  were  true,  we  could  use 
it  to  an  advantage.  Recall  the  modified  post-processor  uses  more  computer 
time  than  the  original.  If  we  could  model  the  Best  Case  results  by  using  the 
3  day  stay  time,  we  could  simulate  the  use  of  the  Attrition  Rate  Model  while 
saving  on  computer  expense. 
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4.3  OTHER  SIMULATION  MODELS 

The  post-processor  presented  in  Appendix  D  is  designed  to  work  only 
in  the  SIDAC  system.  However,  there  are  currently  a  number  of  other  models 
within  the  defense  community  which  are  employed  to  estimate  strategic  fallout 
casualties.  Among  these  are  CIVIC,  COBRA,  READY,  and  RISK  II.  Through 
appropriate  modifications  these  programs  could  be  made  to  handle  the 
Attrition  Rate  equations.  In  this  way  the  Attrition  Rate  Model  of  casualty 
prediction  could  be  more  widely  available  for  use  throughout  the  community. 

Because  modifications  to  the  SIDAC  post-processor  were  straight¬ 
forward  (although  by  no  means  trivial),  there  is  every  reason  to  believe  that 
a  similar  effort  could  be  made  to  modify  the  casualty  prediction  sections  of 
the  other  simulations.  For  example,  in  CIVIC  (Civilian  Vulnerability 
Indicator  Code,  (8))  the  Attrition  Rate  equations  should  be  inserted  in 
Overlay  (5,0),  the  Initial  and  Fallout  Effects  Damage  Assessment  Overlay,  and 
more  specifically,  in  Secondary  Overlay  (5,4),  entitled  EVAL5,  Assess 
Casualties  and  Fatalities  from  Individual  Weapons  -  Initial  and  Fallout 
Effects.  Modifications  to  other  simulation  models  should  be  quite  similar. 

4.4  ALTERNATIVE  DATA  BASES 

The  methodology  outlined  in  Section  2  for  the  construction  of  the 
Attrition  Rate  model  is  not  limited  to  one  specific  data  base.  For  this 
study,  data  were  extracted  from  readily  available  investigations  of  U.S. 
disasters  and  behavior.  Thus,  the  model  is  most  applicable  to  U.S.  popula¬ 
tions. 

If  a  data  base  of  Soviet  disasters  were  available,  one  could 
reperform  the  analysis  to  obtain  a  similar  set  of  equations  based  on  Soviet 
psychological  responses.  This  data  would  more  accurately  reflect  those 
psychological  aspects  which  differ  in  American  and  Soviet  societies.  With 
this  data  base  one  could  feel  more  confident  in  making  statements  about  the 
Soviet  reaction  to  a  nuclear  disaster. 

This  technique  is  not  limited  to  fallout  shelter  studies.  If  an 
appropriate  data  base  exists,  subjects  such  as  industrial  production  under 
adverse  circumstances  or  recovery  after  stress  could  be  investigated. 
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APPENDIX  A 

THE  SIDAC  INSTRUCTION  SET 

Table  A-l  contains  the  set  of  instructions  used  by  CCTC  to  gener¬ 
ate  the  SIDAC  run.  It  is  presented  here  without  comment.  It  may  be  of 
use  to  those  familiar  with  SIDAC  who  desire  to  reproduce  the  results  presen¬ 
ted  in  this  document. 
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Table  A-l.  The  SIDAC  instruction  set  (continued), 
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APPENDIX  B 

CASUALTY  CALCULATIONS 

This  appendix  supplements  the  material  in  Section  3.3.  It 
describes  the  steps  required  to  compute  casualties  in  the  SIDAC  post¬ 
processor. 

The  first  step  is  to  compute  the  proportions  of  the  population 
at  risk  which  become  prompt  fatalities,  prompt  injuries,  fallout  fatali¬ 
ties,  and  fallout  injuries. 

By  assumption,  the  proportion  of  prompt  fatalities  (PFAT) 
equals  the  probability  of  severe  damage  =  1-SPROBS,  where  SPROBS  is 
provided  in  the  Answer  File.  Similarly,  the  proportion  of  prompt  casual¬ 
ties  (PCAS)  ( i . e . ,  fatalities  and  injuries)  equals  the  probability  of 
moderate  damage  =  1-MPROBS.  Therefore,  proportion  of  injuries  =  PCAS-PFAT. 
(SPROBS  and  MPROS  were  previously  defined  in  Section  3.3.1). 

The  fallout  casualty  calculations  are  more  complicated.  For 
each  of  the  7  day,  30  day,  and  180  day  time  periods,  a  "mid-lethal  dose" 
(MLD)  and  a  "mid-casualty  dose"  (MCD)  are  provided.  (Here  casualty  refers 
to  injuries  alone).  Associated  with  these  main  doses  are  their  standard 
deviations,  "standard  deviation  for  lethality"  (SDL)  and  "standard  devia¬ 
tion  for  casualty"  (SOC),  respectively.  For  further  information  on  radia¬ 
tion  doses,  see  (6). 

Recall  that  MAXDOS  is  the  maximum  biological  dose  which  would 
be  received  in  a  given  environment  if  no  protection  existed  (PF=1).  To 
scale  MAXDOS  appropriately,  we  use  that  shielding  value  (SHLVAL)  such 
that  MAXDOS  x  SHLVAL  =  actual  dose  received  by  a  sheltered  individual. 

This  shielding  factor  is  merely  the  reciprocal  of  the  effective  protec¬ 
tion  factor.  For  each  data  record,  there  are  seven  values  of  SHLVAL, 
one  for  each  PF  as  shown  in  Table  3-5.  The  percentage  of  people  corre¬ 
sponding  to  a  given  PF  is  SHLPCT. 

Mathematically,  we  define 

^  DRM(1)/PF1  +  DRM(2)  x  (XMULT( 1 )/PFl  +  XMULT(2)/PF2) 

DRM(l)  +  DRM(2) 


65 


where 


DRM(l) 

=  1.0  -  ITT-0*2 

DRM(2) 

=  ITT-0’2  -  JTT"0,2 

PF1 

=  shelter  protection  factor 

PF2 

=  outside  protection  factor  (assumed  equal 

to 

1/.65) 

XMULT(l) 

=  proportion  of  time  from  ITT  to  JTT  spent 

inside 

shelter,  and 

XMULT(2) 

=  l-XMULT(l)  =  proportion  of  time  from  ITT 

to 

JTT 

spend  outside  shelter 

(See  Section  3.3.2, 

for  discussion  of  ITT  and  JTT.) 

Further  details  on  these  quantities  may  be  found  in  the  Appendix  of  (4). 

As  stated  above  DOSE  =  MAXDOS  x  SHLVAL  is  the  actual  maximum  dose 
received  by  an  individual  in  the  given  fallout  field  and  shelter.  We 
must  now  determine  what  proportion  of  the  population  at  risk  this  dose 
kills  or  injures.  This  is  accomplished  using  the  normal  probability 
distribution. 

Denote  by  $(x)  the  cumulative  normal  probability  at  x.  Let  I 
be  the  index  counting  the  seven  possible  PFs  as  shown  in  Table  3-5.  Then 
the  proportion  of  fatalities  caused  by  the  radiation  dose  DOSE(I)  is 

FAT  =  2  SHLPCT { I )  x  $  ('D0SE,(,1.?  — MLD-) 

1=1  \  SDL  / 


where 


MLD  =  mid-lethal  dose  and  SDL  is  its  standard  deviation. 


Similarly,  the  proportion  of  injuries  caused  by  DOSE(I)  is 


7 

CAS  =  2 

1=1 


SHLPCT(I)  x  $ 


DOSE f I )  -  MCD 
SDC 


\ 


l 


where 


MCD  =  mid-casualty  dose,  and  SDC  is  its  standard  deviations. 

Note  that  MLD,  MCD,  SDL,  and  SDC  are  different  for  each  time 
period  (7  day,  30  day,  or  180  day).  Thus,  this  whole  calculation  is 
repeated  three  times. 

We  now  have  the  four  casual ity  proportions  PFAT,  PCAS,  FAT,  and 
CAS.  The  total  population  at  risk  (capacity  of  the  shelter)  is  CAP. 
Following  steps  3-7  as  outlined  in  Section  3.3.4,  we  calculate  actual 
casualties  as  follows: 

NPF  =  Number  prompt  fatalities  =  PFAT  x  CAP 

NFF  *  Number  fallout  fatalities  =  FAT  x  ( 1 -PFAT )  x  CAP 

NPI  =  Number  prompt  injuries  =  (PCAS-PFAT)  x  (1-FAT)  x  CAP 

NFI  =  Number  fallout  injuries  =  (1-PCAS)  x  (1-FAT)  x  CAS  x  CAP 

NHP  =  Number  healthy  population  =  CAP-(NPF+NFF+NPI+NFI ) 

This  completes  the  computation  for  one  data  record.  We  repeat 
this  entire  procedure  for  each  data  record  and  sum  the  results. 

To  illustrate  this  procedure,  consider  the  following  example. 

Let  MPROBS  =  .5  and  SPROBS  =  .7.  This  means  that  probability  of  moderate 

damage  =  PCAS  =  .5,  and  probability  of  severe  damage  =  PFAT  =  .3.  Let 

CAP  =  1000,  MAXDOS  =  4000.  Assume  that  for  the  time  period  in  question, 

MLD  =  500  MCD  =  200 

SDL  =  200  SDC  =  50. 
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For  simplicity,  assume  the  following  shelter  distribution: 


PF 

PERCENT 

200 

20 

50 

50 

10 

30 

We  compute  casualties  resulting  from  a  3  day  stay  time 
by  18  days  of  2/3  sheltering. 

First,  we  compute  three  values  of  SHLVAL,  one  for  each 
ITT  =  end  of  period  of  100%  sheltering  =  3  days  =  72  hours.  JTT 
partial  sheltering  =  21  days  =  504  hours. 


DRM(l) 

= 

1  -  ITT'0'2  =  .575 

DRM(2) 

= 

ITT-0’2  -  JTT-0*2  =  .137 

PF1 

* 

200,  50,  and  10,  respectively 

PF1 

= 

1/.65  *  1.538 

XMVLT(l) 

= 

2/3  =  .667 

XMVLT(2) 

= 

1/3  =  .333 

Therefore,  for  PF1  =  200, 

SHLVAL (1 )  -  .575/200  +  ,137( .667/200  +  .333/1.538) 

.575  +  .137 

=  .046 

Similarly,  for  PF1  =  50,  SHLVAL(2)  =  .058,  and  for  PF1  =  10, 
SHLVAL(3)  =  .135. 


followed 

PF. 

=  end  of 
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We  now  compute  FAT  and  CAS  by  filling  in  Table  B-l.  By  summing 
the  two  indicated  columns,  we  find  FAT  =  .23  and  CAS  =  .74,  i.e.,  23%  of 
the  at  risk  population  are  fatalities  and  74%  are  casual ities.  From 
before,  we  had  PFAT  =  .30  and  PCAS  =  .50.  We  now  compute 

NPF  =  .30  x  1000  =  300  prompt  fatalities 

NFF  =  .23  x  .70x1000  =  161  fallout  fatalities 

NPI  =  (.50  -  .30)  x  (1  -  .23)  x  1000  =  154  prompt  injuries 

NFI  =  (1  -  .50)  x  (1  -  .23)  x  .74  x  1000  =  285  fallout  injuries 

NNHP  =  1000  -  (300  +  161  +  154  +  285)  =  100  uninjured  people. 

This  data  record  has  been  completed;  we  would  now  get  another 

record,  compute  the  same  quantities,  and  add  them  to  these  results. 


69 


Table  B-l .  Fallout  casualty  calculations. 
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APPENDIX  C 

THE  ORIGINAL  SIDAC  POST-PROCESSOR 
Table  C-l  contains  a  listing  of  the  original  SIDAC  post-processor 
as  written  by  the  Command  and  Control  Technical  Center.  What  follows,  while 
not  a  complete  documentation,  is  intended  to  serve  as  a  guide  to  the  program 
segments. 


LINE  NUMBERS  IN  TABLE  C-l 


1 

-  30 

31 

-  72 

73 

-103 

105 

-119 

120 

-134 

135 

-170 

171 

-238 

239 

-342(end) 

FUNCTION 
Initial ization 

Initial  data  acquisition,  continued 
initial ization 

Read  a  record  and  assign  values  to 
variables 

Default  values  for  PFs  and  percen¬ 
tages  —  urban  population 
Default  values  for  PFs  and  percen¬ 
tages  —  rural  population 
Additional  initialization  and  error 
calls 

The  bulk  of  the  computations  --  see 
further  comments  below 
Output  and  bookkeeping  routines. 


As  noted  above,  lines  171-238  comprise  the  bulk  of  the  computations. 
To  aid  in  comprehension,  some  of  the  key  variable  one  identified  here.  The 
actual  flow  of  the  computations  is  fairly  clear. 

There  are  many  do-loops  which  run  from  1  to  7.  (See  lines  171, 

184,  187,  190,  etc.)  .  These  index  the  seven  different  PFs  per  environment. 
Do-loops  from  1  to  3  index  the  three  different  assessment  times,  7  days,  30 
days,  and  180  days.  PR1  is  the  inside  protection  factor;  PR2  is  the  outside 


i 

{.  m**  -  "  ■ 
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protection  factor;  set  equal  to  1/.65.  The  proportion  of  partial  shelter¬ 
ing  spent  indoors,  denoted  XMULT(l),  is  equal  to  1-XMULT  (2).  ITT  is  time 
at  end  of  complete  sheltering;  JTT  is  time  at  end  of  partial  sheltering  and 
marks  beginning  of  evacuation.  DRM(l)  and  DRM(2)  are  dose  rate  multipliers 
for  each  time  period.  SHLVAl  is  the  reciprocal  of  the  effective  protection 
factor. 

SMLD  (I),  SIGL(I),  SMLD(I),  and  SIGC(I)  are  mid-lethal  dose,  its  stan¬ 
dard  deviation,  the  mid-casualty  dose,  and  its  standard  deviation,  respec¬ 
tively,  for  the  assessment  time  indexed  by  I  *1,  2,  or  3.  The  actual 
values  they  assume  (in  lines  224-231  and  lines  15-18)  are  from  reference 
(6).  (MAXDOS)  x  (SHLVAL)  is  the  actual  dose  received.  CUMN  is  a  sub¬ 
routine  which  computes  the  cumulative  normal  distribution.  FAT  and  CAS 
store  the  percentage  of  fatalities  and  casualties,  respectively.  In  lines 
240-246,  these  percentages  are  converted  to  actual  population  counts. 

To  run  this  program,  the  Answer  File  must  be  available  as  device 
number  25.  Device  5  is  the  card  reader,  so  all  reads  to  device  5  must  find 
the  data  located  after  the  source  code. 


Table  C-l.  The  original  SIDAC  post-processor 


1 

2 

3 

4 

s 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

6020/ 78W 


91  MENS  ION  I  PR Of (2,5,7) , I SHPCT (5 , 7 ) , I NUL T ( 2, 5 > , I SHL V ( 7) , 
11SHLP(7),DRM(2),ITX(4),IPCTF(5> 

DIMENSION  XMULT(2,5) 

CHARACTER  I  CC  *2 ( 5) 

CHARACTER  TITLE1*20/TIT|.E2*20/TITLE3»20/TITLE4*S 
CHARACTER  JNAM«5,T0P0S*8,SPAC1»4»SPAC2*3 
CHARACTER'  I  SUB* 1 ,S  VSU3*  1 

DIMENSION  AL1N(3,6,3),PCT(7),ATEN(6,3),IFAT(3>,ICAS(3),FAT(3), 

1  CAS(3),C0M(3S) 

CHARACTER  IRG»1  , SVRG*1 , 1 TI HE* 3 < 3) , I N AM*5 (2 ) , 1 C AT  *  1 

CHARACTER  ICTY*2/*  ,/,SVeTY*2/*  */ 

DATA  I T X  724,72,240/723/ 

DATA  I TI ME /  1  7*,*  30*, *180'/ 

DIMENSION  SHLVAL(16)/SHLACT(16)/SMLD(3)/SNCD(3),SIGL(3)/SIGC(3) 
DATA  SMLD/1000.,0.,453. / 

DATA  SMCD/200.,215.,433./ 

DATA  SIGL/200.,1 .,1 35./ 

DATA  SIGC/60.,64.5,123./ 

DATA  SHLVAL/.01,.075,.15,.2,.3,.55,.7, 

(.1 ,.15, .2, .25, .3, .4, .55, .6, .7/ 

DATA  SMLPCT/.05,.12,.045,.385,.09,.21 ,.1, 

1.307,. 00 3 5,. 036,. 01 05,. 06 5,. 075,. 25 3,. 25,. 3/ 

DIMENSION  XSHLV(8)  ,KSH(.P(400> 

CHARACTER  KCTY*2 (50) 

DATA  XSHLV/200. ,130.,53.,20.,10.,5.,3.,1 .4/ 

DATA  ISHUV/200,100, 50,20, 10,5,3/ 

REAL  MPROBS 
IPAC*0 
REWIND  25 
REWIND  10 
1*1 

KCTY(I)*"  " 

1011  REAO(10,1012,EN0*101S)<CTY(I),(KSHLP(8»(I-1>*J),.l*1,8) 

1012  FORMAT <A2, 812) 

1*1*1 

GO  TO  1011 

1013  CONTINUE 

READ  (5,100)  TITLE1/TITLE2, TITLES, TITLE! 

100  FORMAT  (3A20/A3) 

READ  (5,449,EN0*21 21)  I TT , JTT 
449  FORMAT  (14,14) 

IF(JTT.EO.O)JTT*504 
IF (JTT.SE.ITT)GO  TO  3131 
KTT*J  TT 
J  T  T*  t  T  T 
I T  T *<  T  T 
3131  CONTINUE 

DO  46  J  *  1 , 5 

READ  (5/30»€N0*46)  I C C ( i ) , I HULT( 1 , J > , I  MULT ( 2 , i ) , 

1  (IPR0F(1,J,<),K*1,7),(IPR0F(2,J,K),X*1,7), 

1(ISH*CT(J,K),K*1,7),I°CTF(J) 

30  FORMAT (  A2, 213, 1413, 712,  ID 

46  CONTINUE 
>121  CONTINUE 
I SUB»1H 
S  N  $  U3  *  1  H 
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Table  C-l.  The  original  SIDAC  post-processor  (continued) 


57 

NU*1  000000. 

58 

JNAN, 

•urban* 

59 

S  9 A  Cl *4  H 

60 

SPAC2*3H 

61 

SVCTY*2H 

62 

SVRG«1H 

63 

TOPOS»8H 

64 

I EO  F*0 

65 

ILIN* 

60 

66 

I  F  L  A  3 

•0 

67 

00  21 

11*1,3 

68 

00  21 

12*1 ,6 

69 

00  21 

13*1,3 

70 

AIIN< 11/12, I3>«0. 

71 

A  T  E  N( 12, 13) *3. 

72 

21 

CONTINUE 

73 

6 

HEAD 

( 2  5 , E  N0*9  8 ) 

CON 

74 

NPR03S*C0N (4) 

75 

SPR03S*C0N(5> 

76 

CALL 

BYTE! COR, 163 

1CAP,1,6) 

77 

CALL 

8YTE(C0N,183 

I RG, 1 , 1 ) 

78 

CALL 

BYTE (CON, 1 81 

ICTY, 1,2) 

79 

CALL 

BYTE(C0N,1 89 

I  SUB , 1,1) 

80 

CALL 

BYTE(C0N,1 71 

IFAT(1),5,2) 

81 

CALL 

BYTE (CON, 175 

IFAT(2),5,2> 

82 

CALL 

BYTE(C0N,1 79 

1 F  AT ( 3) ,5,2) 

83 

CALL 

BYTE ( C0N,1 69 

ICAS(1),5,2) 

84 

CALL 

BYTE ( CON  ,  1 73 

I  CAS ( 2) ,5,2) 

85 

CALL 

B  Y  T  E ( CON ,177 

I C  AS ( 3) ,5,2) 

86 

CALL 

BYTE ( CON, 1 84 

I  CAT, 1 , 5) 

87 

CALL 

8YTE(CON,043 

NAXOOS, 1 ,6) 

88 

CALL 

BYTE ( CON ,037 

ID NODS/1, 6) 

89 

CALL 

BYTE(C0M,193 

ISHLI/(1  >,4,3> 

90 

CALL 

BYTE (CON, 1 96 

ISHLV(2),4,3) 

91 

CALL 

BYTE  <  C0N,1 99 

I  SHL  V  ( 3  >  ,4,  3) 

92 

CALL 

8YTE (CON, 202 

ISHLV(4),4,3) 

93 

CALL 

BYTE (CON, 205 

ISHLV(5),4,3) 

94 

CALL 

BYTE (CON, 208 

ISHLV(6),4,3> 

95 

CALL 

BYTE(CON,223 

ISHLV(7>,4,3) 

96 

CALL 

BYTE ( CON, 1 91 

ISHLP(1 ),5,2) 

97 

CALL 

BYTE(C0N,211 

ISHL»(2),5,2> 

98 

CALL 

BY  TE ( CON, 2  1  3 

ISHL»(3),5,2) 

99 

CALL 

BYTE( CON, 21  5 

ISHLP(4),5,2> 

1  00 

call 

BYTE (CON, 21 7 

I  SHL" (5 ) ,5,2) 

1  01 

CALL 

8YTE(C0N,219 

ISHL°(65,5,2) 

1  02 

CALL 

BYTE(C0N,221 

ISHLP(7),5,2) 

1  03 

IF  (  ICAT.EQ  ."75399 

".OR. ICAT.EQ.” 

1  04 

I  SHLV ( 1 >  *800 

1  05 

ISHLY(2>*250 

1  06 

ISHLt/(3)*150 

1  07 

ISHLrf(4)»100 

1  08 

l$HLi/(5)*40 

1  09 

I  SHL9  <6 ) *20 

1  10 

I SHLV ( 7) *1 0 

1  1 1 

1SHL9  <1 1*1 4 

1  12 

I  SHL9 (2  >«7 

750990 ”>GO  TO  251 


Table  C-l.  The  original  SIDAC  post-processor  (continued) 


13 

14 

15 

16 

17 

18 

19 

20 


22 

23 

24 

25 

26 
27 


1  29 
1  30 
1  31 
I  32 
1  33 
I  34 
i  35 
1  36 
1  37 
1  38 
1  39 


1  42 
I  43 
I  44 
1  45 
1  46 
1  47 
I  48 
I  49 
I  50 
151 
I  52 
I  53 
I  54 
I  55 
I  56 
I  57 
I  58 
I  59 
I  60 

61 

62 

63 

64 

65 

66 

67 

68 


I$HLI»<3)«7 
ISHLP<4)»14 
I SHUP  <5 ) *28 
ISHL®<6>«1 7 
ISHLP<7)*1 3 
SO  TO  252 

251  ISHLV(1)«40 
I S  HU V  <21*35 
I  SHUV <  3) »30 
I SHLV (4) »2  5 
I SHLV ( 5) ■?  0 
I$HLV<6)«1 5 
I S  HLV ( 7 ) *  1 0 
ISHL»<1 >«2 

I SHUP ( 2 ) *  5 
ISHLP(3)*1 5 
1 SHL°  (  4  )  *8 
ISHU#<5>*22 
I  SHU* 16) *6 
ISHLP<7)«42 

252  CONTINUE 

if<  icTr.ea.,*us")S3  to  ioi? 

DO  1314  JJ«1,50 

If <KCTV(4J).6a."  ">S0  TO  1017 

1314  IF<ICTY.Ea.KCTY<JJ)>SO  TO  1015 

1015  00  1316  KK«1#7 
ISHLV<KK)«XS4Ly  <)U> 

1016  ISHL»<KK)«KSHUP< ( i i -1 ) » 8» K X) 

1017  CONTINUE 

IF  (I  FLAG. EO. 1 )  GO  TO  52 
SVCTY*ICTY 
SVRG* I  PS 
SVSU3* I $U8 
52  CONTINUE 

CAP*FUOAT{  I C A P ) 

IF(ICAT.Ea."75100".OR. ICAT.E8."751000’')CAP*CAP«1030. 
IF(IC4T.EQ."75399".0«. I C A T . EO . " 750990" ) C AP»C AP *  1 00 . 
DO  23  J  *1 »  3 
F4TU)«IFA  TCJ  > • . 01 
CAS  (J  >*ICA$0  >*.01 
23  CONTINUE 
I  N  T  *1 
I  N  F  «7 

IF  <I»AC.E8.5)  I  NT*8 
IF  <IPAC.E8.5>  I NF ■ 1 6 
00  191  13*1.3 
CAS (I 3>*0. 

FAT(I3)*0. 

191  CONTINUE 
11*1 

44  i f< n .ea.6)  go  to  201 

I F< ICCC 1 1) .E8.ICTY.0R. ICC (II >  .EQ. 'XX')  SO  TO  146 

I  I  *  I  I  A  1 
GO  TO  44 

201  PAINT  202 

202  FORNAT  <15X.  mill,  'CARO  XX  HISSING') 


Table  C-l.  The  original  SIDAC  post-processor  (continued). 


1  69 

CALL  EXIT 

1  70 

146 

00  143  LL*1,7 

1  71 

143 

SHLPCT (LL)*FLOAT (ISHLP(LL) >/100. 

1  72 

C 

THE  CALCULATION  USES  ONLY  THE  PCT  61ST.1N  THE  DATA 

RECORD. HOWEVER 

1  73 

C 

ISHPCT  HAS  SEEN  READ  IN  AND  PGM  CAN  BE  MODIFIED 

1  74 

c  •** 

IT  IS  ASSUMED  THAT  THERE  WILL  BE  ONLY  ONE  PCT  DIST. 

FOR 

BOTH 

1  75 

c  «•* 

(POSSIBLE)  SETS  Of  PROTECTION  FACTORS (IPROF) . 

1  76 

c  »** 

THE  FORM  OF  THE  CALCULATION  ASSUMED  FOR  THE  SECOND 

TIME 

PERIOD 

1  77 

c  **» 

IS  1ST  NUL  T* ( 1 » / 1ST  PROT.  FACT.)  ♦  2ND  HULT*(1./2ND 

P.F 

.). 

1  78 

C  »'* 

ALSO  IT  IS  ASSUMED  THAT  1ST  PERIOD  USES  THE  PROT.  FACT. 

IN  THE 

1  79 

c  *•* 

DATA  BASE  RECORO. 

1  80 

c  *»* 

1  81 

c  *•* 

I • C  TF  IS  THE  FLAG  TO  USE  NEW  PCT  DIST.  OR  NOT. 

1  82 

IF(I®CTFCI  D.EO.0)  GO  TO  144 

1  83 

00  145  LL*  1 , ? 

1  84 

145 

SHLPCT(LL)*fLOAT<ISHPCT(lI,LL))7iaO. 

1  85 

144 

IFdI.LT. 6)  GO  TO  147 

1  86 

DO  148  J  J  ■  1 , 7 

1  87 

148 

SHLVAL<JJ)®F.OAT(ISHL/(JJ))/10D. 

1  88 

GO  TO  991 

1  89 

147 

DO  149  K*1,7 

1  90 

PR  1  "FLOAT  (  IPROF  O,  I  l,K>  ) 

1  91 

PR2*FL0AT ( I PR OF (2, I I ,  <  > ) 

1  92 

IF(IPROF(1,lI,lO.LE.O>  PR1»FL0AT<  ISHLVU)  ) 

1  93 

IFd*ROF(2,II,K).LE.0)  PR2»FL0AT(ISHLV(K)  ) 

1  94 

c 

IF  THE  PROTECTION  FACTOR  IN  THE  DATA  BASE  RECORD  AND  THEPROT. 

1  95 

c  ••• 

FACT.  IN  THE  CHAN6E  CARO  ARE  BOTH  0,  IT  IS  ASSUMED 

THAT 

THE 

1  96 

c 

PROT.  FACT.  •  1./.65 

1  97 

1F(I*R0F<1,1I,IC).IE.D.4ND.ISHLV<X).LE.0>  PR  1*1 .7.65 

1  98 

IF<  I*R0F(2,II,K) .LE.O.ANO.ISHLVOO  .LE.O)  »R2*1  .7.65 

1  99 

IF  dPROFd  ,11  ,K)  ,E8.  1)  PR1«1./.65 

2  00 

IF (  I*ROf (2,11 ,K> .:8.  1)  PR2*1 . 7 . 65 

201 

IF(DRM(1)®ORM<2).LE.0)  SHLVAL  CX>*.65 

202 

X  MULT (1,11 )*FLOAT( IMULT (1  ,1 1) )/100. 

203 

XMULT (2, II )*FLOAT< IMULT (2,1 1 )>7100. 

204 

IF<  I  MULT (1 ,1 r 1.E 0.0. A  NO. IMULT (2,1  I) .E0.3)  XMUL T ( 1 , 1 1 > * 1 

.  0 

205 

DRM (1 ) *1 -I TT*«<-.2) 

206 

ORM(2)*lTT«*(-.2)-JTT*»(-.2) 

207 

IF((ICTY.E9."JA".0R.ICTY.E8."TW',).AND,K.EB.1)PR1*1. 

4 

2  08 

IF<ICTY.E0."TW".AN0.X.ER.1)SHLPCT(1)*.3 

209 

IF  <IC  TY.EQ."J A".AN0.X.E8. 1 ) SHLPCT (1 )*.4 

2  10 

FAC1*0RM(1  )*<1./PR1) 

2  11 

FAC2*XMULT (1,II)*(1./PR1) 

2  12 

FAC  3* XMULT (2, II >*(1./PR2> 

2  13 

FAC4*0RM(2)  •(FA:2»FAC3) 

2  14 

FAC5*0RM(1 )>DRM<2) 

2  15 

F  A  C6*  F  AC  1 ♦ FAC  4 

2  16 

SHLVAL0C)*FAC67FAC5 

217 

149 

CONTI NUE 

2  18 

991 

CONTINUE 

2  19 

I  NT*1 

2  20 

I  N  F  *  7 

2  21 

DO  999  11*1,3 

2  22 

00  999  I»*INT,INF 

2  23 

SMLD(2)*1000. 

2  24 

SIGL(2>*200. 

76 


Table  C-l.  The  original  SIDAC  post-processor  (continued) 


225 
2  26 
22? 
223 
2  29 
2  30 

231 

232 
2  33 
2  34 
235 
2  36 
237 
2  38 

239 

240 

241 

242 

243 

244 

245 

246 

247 
2  48 
2  49 
2  50 

251 

252 

253 

254 

255 
2  56 
2  57 

258 

259 
2  60 
2  61 
2  62 

263 

264 

265 
2  66 
2  6? 
2  68 
2  69 
2  70 
2  71 
2  72 
2  73 
2?4 
2  75 
2  76 
2  77 
2  78 
2  79 
2  80 


IF(M4XOOS»SHLVAL(I2>  .LE.400.)SNLD(2)»540. 

IF(NAXDOS*SHU/AL (I 2 ) ..€ .4 00 . ) SI GL <2 ) «1 62 . 

IF  <M4X00S<SHIVAL<1  2)  .  G T . 400 . . AND. MAX DOS  * S HU/ AL  <I2)  .16.1  300.) 

%  SNLD(2)  »?20. 

IF  (MAX00S*SMLVAL (12) .GT.4  00..AND.NAXDOS*$HLt/AL (I  2) .LE. 1  300. > 
t  SIGL(2)»3?0. 

FAT<I1)aFAT(I1)ACUMN((NAX00$»SHLVAL(I2)-$MLD(I1))/SIGL(I1)> 

t*SHLPCT(I2) 

CASU1>«CAS(I1>*CUNN(<NAXDOS<»SHLVAL<I2)-$NCO<I1))/SIGC<I1>> 

t*SHL»CT(I2) 

I  F  (  (NAXDO$*SHLVALC 12)> .  £8.0)  FAT{  II  )»0. 
IF((NAXD0S«SHU/AL<I2)).E8.0>  CAS(11)»0. 

999  CONTINUE 

DO  22  I  »1  #  3 
AT£N-(6#I)*CA° 

ATEN(5»I)«(1.0-SPR09S>»CAP 

ATEM<1/I)»SPR08S*CAP*F4T(I) 

ATEM(2/1)*<SPR09S-NPR03S)«<1.0-FAT(I))*CAP 

ATEM(3#I )«HPR08S*(1 . 0-F  AT  < I > ) *C AS ( I > » C AP 

ATEN(4#t )»CAP-(ATEM(5,I )AATEM(1 /I )*ATEM(2#I >44  TEN (3/ I ) > 

IF  <ATEN(4#I) .L  T.3. )  ATEN(4,l>»0. 

22  CONTINUE 

IF  (I  FLAG. E8. 0)  GO  TO  4 

IF  (ICTY.NE.SVCTY.0R.IRG.NE.SVRG.0R.ISU8.ME.SVSUB)  Il«1 
IF  (ICTY.NE.SVCTY.OR.IRG.NE.SVRG.OR.ISUB.NE.SVSUB)  GO  TO  7 

4  IFLA3«1 

DO  5  Jl»1 #3 
00  5  J2»1,6 
DO  5  J 3»1 #  3 

ALINI  J1#J2,J3)»AUN(J1  #  J2,J3>4ATEH(J2,J3) 

5  CONTINUE 

IF  (IE0F.EQ.1 )  GO  TO  7 
GO  TO  6 

7  1LIN«IUN*4 

IF  (IUN.LE.53)  GO  TO  3 
33  PRINT  69 

69  FORMAT  ( 1 H  1  ) 

PRINT  61 
PRINT  62/TITLEl 
PRINT  63 / T  ITLE2 / JNAN 
°R  I  NT  64/  TITLES 
PRINT  65/TITLE4 
PRINT  76 
PRINT  68 
PRINT  66 

61  FORMAT  <  TZ / * • •  RED  ON  BL U E * / 40X / 1 POPUL AT  1  ON  ASSESSMENT*) 

62  FORMAT  <  T2  /  * • *  C  AS  E  /  SC  E  NA  R  10  :  */A20) 

63  FORMAT  CT2/*#*  SPEC  INSTR:  * / A20, 26 X / A5 ) 

64  FORMAT  <T2#*»»  SPEC  INSTR:  * / A20/ 26X / 5 < 1 H- ) ) 

65  FORMAT  (T2/*»»*#1X/A3#1 X/ 'ASSESSMENT* ) 

68  FORMAT  ( T37/' FATAL  I  TIES */ 36X/ *1 NJURIE S * ) 

66  FORMAT<  T22,43< 1M-) /2X/43 ( lH-)/6X, ’RESIDUAL '/ 3X, *  ASSESS  * /// 

3  1 2  X/ ' TO  T  P0P*/6X/2(*PR0MPT*/8X,* FALLOUT* ,1DX, ’TOTAL* /9X), 

t  2X/* P0P*/7x, *TI ME './»12x #7(1 H-).2X #7(13(1 H-)/2x)/6(1H-)/// 

S  1X/'REG*/2X/*CTRY*/4X/*MIL'#7X/6(*MIL'/3X/'PCT*/6X)» 
i  ’MIL*»3X#'PCT*,4X,*DAYS*#/#1X#3(1H-)»2X#4(1N-),3X#5(1H->,3X, 
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Table  C-l 


The  original  SIDAC  post-processor  (continued) 


281 

3  7<8<lH-)/1X/4(lH-)/2x>/6(1H->///> 

282 

I  L I N*  1  7 

2  83 

8 

03  13  K2*1/6 

2  84 

03  13  K3*1/3 

28$ 

ALIN.(I1/K2/K3)*ALIN<I1,K2/K3>/NIL 

2  86 

1  3 

continue 

287 

GO  TO  <40/41/421/11 

2  88 

43 

ENCOO£(TOPOS/50)  SVRG/SPACI/SVCTY 

2  89 

50 

FORMAT  (A1/A4/A2/*  *) 

2  90 

GO  T3  75 

2  91 

41 

ENCODE  (TOPOS/51)  SVR3.SPAC2 

2  92 

51 

FORMAT  (Al/A3»*ALL  *) 

2  93 

GO  T3  75 

2  94 

42 

ENCODE  (TOPOS/57) 

295 

57 

FORMAT  ( 'WW*/2X/ ' ALL  *  ) 

2  96 

75 

DO  9  13*1/3 

2  97 

T*A*ALIN(I1/5/I3)*ALIN(I1/1/I3) 

298 

TIN*AH*.(I1/2/l3)*ALIN(I1/3/I3) 

2  99 

00  10  K*1  /  5 

3  00 

PCT  (<)*(AL  IN<  11  /*/  13)/  ALINU1  /6/I3)  )*100. 

301 

10 

CONTINUE 

302 

®CT(6)*<TFA/ALIN(I1/6/I3))»100. 

303 

PCT (7)* (TIN/ALINU 1/6/13) >*100. 

3  04 

PRINT  1/  TOPOS  . ALI N< I 1/6/13) /ALINd 1 /5/I 3 > *PCT ( 5) / 

305 

3  ALIN(I1/1/I3)/PCT(1)/TFA/PCT(6)/ALlN(I1/2/I3>/*CT(2)/ 

306 

3  ALIN(I1/3/I3)/PCT(3)/TIN/PCT(7)/ALIN(I1/4/I3)*»CT(4)/ITIME(I3) 

307 

1 

format  (2X/A5/  1X/F5.3/2X/7(F8.3/1X,F4.1,2x>,1X,A3> 

308 

9 

continue 

309 

IF  (I1.E0.2)  PRINT  76 

310 

IF  (I1.EQ.2)  I L I N* I L I N4  1 

311 

PRINT  76 

312 

76 

FORMAT  </) 

313 

00  11  12*1/6 

314 

DO  11  13*1/3 

315 

AlIN< I1/I2/I3)*0.0 

316 

1  1 

CONTI NUE 

317 

IF  (I1.E0.1)  SVCTY*ICTY 

318 

IF  (I 1 .E9. 3. AND. IEOF. ED. 1 )  GO  TO  99 

3  19 

IF  (I1.EQ.2)  SVRG»IRG 

320 

IF  (I1.E0.2)  GO  TO  12 

321 

IF  (I1.E0.3)  SVSU3* I SU3 

322 

IF  (I1.E0.3)  I L I N*  60 

323 

IF  <1 1 ,E9. 3)  J N AM  * ‘RURAL  * 

3  24 

IF  (I1.EQ.3)  GO  TO  4 

325 

IF  <1 RG.NE .SVRG.OR. I SU9.NE. SVSUB. OR.  IEOF  .£0.  1  >  11*2 

326 

IF  (IRG. NE.SVR3. OR. ISJ9.NE. SVSUB. OR. IEOF. EB. 1 )  30  TO  7 

327 

GO  TO  4 

328 

1  2 

IF  (ISU8.NE. SVSUB. OR. IEOF. E3.1)  11*3 

329 

IF  (I SU9.NE. SVSUB. OR . IEOF ,EQ. 1 )  GO  TO  7 

3  30 

GO  T  3  4 

3  31 

99 

I  E  0  F  *  1 

3  32 

11*1 

333 

GO  TO  7 

3  34 

99 

STOP 

3  35 

END 
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APPENDIX  D 

THE  MODIFIED  SIDAC  POST- PROCESSOR 


Table  0-1  presents  the  SIDAC  post-processor  as  modified  to 
handle  the  Attrition  Rate  equations.  It  is  slightly  longer  than  the 
original  (361  lines  to  342  lines).  Other  than  some  overall  modifications 
and  simplifications  (which  could  be  equally  well  applied  to  the  orginal) 
the  major  changes  occur  in  the  prime  computational  section,  lines  171  to 
252,  and  at  line  32.  (Other  changes,  such  as  various  initializations,  will 
become  obvious  upon  comparing  the  two  sets  of  code). 

Line  32  reads  the  two  coefficients  from  the  attrition  rate 
equation,  y  =  AA  x  t  +  88.  On  lines  171  to  252,  the  primary  modification 
is  the  addition  of  two  loops  (line  206  and  238)  indexed  from  1  to  N1 ,  where 
N1  is  the  number  of  time  periods  (including  the  "Zeroth")  until  the  entire 
population  has  left  the  shelter.  All  fatalities  and  casualties  are  cal¬ 
culated  separately  for  each  time  period,  and  a  separate  SHLVAL,  denoted 
SHVALX,  is  computed  for  each.  The  dose  received  for  each  group  is  held  in 
STOR.  FATX  and  CASX  contain  separate  casualty  percentages  for  each  group. 
These  separate  figures  are  finally  recombined  in  lines  254-263  by  multi¬ 
plying  by  the  proportion  of  people  in  each  group.  The  remainder  of  the 
program  is  identical  to  the  orginial. 

To  run  this  program,  follow  the  same  procedure  as  in  the  original 
post-processor,  except  a  card  giving  the  values  of  AA  and  BB  must  be 
included  to  be  read  at  line  32. 


Table  D-l.  The  modified  SIDAC  post-processor 


1 

2 

4 

5 
t 
7 
b 
9 


1  0 

1 1 

1  2 

1  3 

1  4 

1  5 

1  0 

1  7 

13 

1  9 

20 

21 

22 

23 

24 

25 

2  o 

27 

2a 

29 

30 

31 

32 

33 

7  7o 

34 

35 

; 

36 

101  1 

37 

1012 

3o 

39 

40 

1013 

41 

42 

100 

43 

44 

449- 

45 

4  o 

47 

43 

49 

50 

SI  31 

5  1 

52 

53 

54 

55 

30 

5a 

46 

DIMENSION  I  PROF ( 2,5,  7) ,  i SHPC  T( 5.7) i I MULT (2/ 5 ) * ISHLVC  7) , 

KISH UP (7), DR  1(2), ITX(4),IPCTF(5) 

DIMENSION  Xi.ULT(2,5> 

CHARACTfcR  I  C  C  *  2  (  5  ) 

CHARACTER  TITL£1*2G,TITLE2  +  2G,TlTLc3+20,TITLt.4  +  3 
CHARACTtR  J  .JA.ri »  5  ,  T  C.PO  S  »  „,  SP  A  C  1  +4,SP4C2*3 
CHARACTER  I  3  Ud  *  1  ,  S  V  j  u  U  "  1 

01 KENS  I  ON  ALIN(3,6,3)/PCT(7),ATEI'(6/3)/IfAT(3)/ICAS(3)/FAT(3)/ 

1  CAS<3)/COM(3o) 

CHARACTER  IhU*1  , jtfHb+1 ,ITIM£*i< 3) /INAM+5 (2) ,IC4T»1 

ChARACTER  ICTY *  2  /  1  ,/,SVCTY*2/'  '  / 

j  A  T  A  I  T  x  /24,72/24u,?2L/ 

DATA  I  T  I  ME/  •  7  1  ,  *  jC'/MoU'/ 

DIMENSION  SHLVAL(1o>,ShlPCT(16),SiiL0(3),S''1CD(3),SIGL(3)/SIGC(3) 
lata  SMLD/1 000., P.,450./ 

DATA  SiiCu/2  0u.,215./4u0./ 

OAT  A  SIGL/200.,1  .  /  1  3  5  .  / 

DATA  SluC/SU. ,04.5/120./ 

D  A  T  A  6  H  L  i  A  w  /  .01  /  .  0  7  a  ,  .  1  3/.2/.3/«^S,  .7  , 
&.1/.15/.2/.25/.3/.4/.55/.6/.7/ 

OATA  ShLPCT/.C5/.12s.u<»S/«3oS/.C9,.21/.1/ 

2  .0C7,.uC55/.G36,.31C5/.l.'o5,.u75/.2  53/.25,  .3/ 

DIMENSION  X  S  H  L  V (  6 ) ,  A  S  H  L  P  (  4  OC ) 

Character  kcty«2(5C) 

DATA  XSHLV/2uJ./1CC.,50./2u./1'J./:'./3./1  *4/ 

DATA  I  She V / 2UO , 1 00 , 50 , 2 0, 1 j , 5 , 5 / 

REAL  M  P  H  J  b  S 
IPAC=0 
REMIND  23 
REwlND  1J 

DIMENSION  SHVALX(16,5G),FATa(3,:>(J>,CASX(5,5j) 

R  E  A  li  (  4  7 , 7  7o  )  *Oo  b 
f  D  R  K  A  T  <  2  F  1  0  .  G  ) 

1=1 

nCTY(I)«" 

R£AD( 10/101 2,tND=1G13)ACTY (I ),(KSHLP(6*( 1-1  )+d),J=1/2) 

FORMAT  <  A  2  /  o  I  2) 

1  =  1  +  1 

GO  TO  1011 
CONTINUE 

READ  (5/100)  T I TLE 1 / T I TL22/TI TLti/ T  ITLt4 
FORMAT  (3A20/A3) 

READ  <5/449, £10=21 21 >  ITT/JTT 
FORMAT  (14,14) 

IF(JTT.EG.U)JTT=534 

IF(jTT.bE.ITT)GO  TO  3  1  3  1 

A  T  T  =  J  T  T 

J  T  T  =  I  T  T 

I  T  T  =  K  T  T 

CONTINUE 

DO  4S  J=1 ,5 

READ  ( 5 , 30, END  =  4o)  I CC ( J ) , IMULT(1 ,J > , IMULT (2/J) , 

&  (IPKUF(1,J,K>,K=1,7),(IPR0F(2,J,K),K=1,7), 

6(ISHPCT(J/A)/A=1/7)/IPCrF(J) 

FORMAT (  A2, 213/1413/712,11) 

CONTINUE 


6020/78  W 
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Table  D-l.  The  modified  SIDAC  post-processor  (continued) 


57 

2121 

CONTINUE 

58 

I  SUd 

■  1  H 

59 

SVSUb* 1 H 

60 

MIL»1 OOCOOO. 

61 

J  NAM* 1  URBAN  * 

62 

SPAC1»4H 

63 

SPAC2*3H 

64 

S  V  C  T  Y  *  2  H 

65 

SVRG*lH 

66 

T  OPQ  S  a  8H 

67 

IcOf *0 

66 

1L1N 

*6U 

69 

IfLAGaO 

70 

D  0  21  11=1,3 

71 

DO  21  12=1,6 

72 

DO  21  13*1,3 

73 

ALI.\i(I1,I2,I3>*0. 

74 

ATEMt 12,13) *0. 

75 

21 

CONTINUE 

76 

6 

READ 

( 2  5, END=9  8)  CON 

77 

MPR0u$*C0M(4) 

78 

SPROBS=COM<5> 

79 

CALL 

UYTE(C0K,1o3,ICAP,1,o) 

60 

CALL 

BYTE  (C0M,1  33,186,1,1) 

81 

CALL 

8YT£(C0N,1a1, ICTY, 1,2) 

62 

CALL 

6YT6<C0«,169,1SU0,1,1> 

83 

CAUL 

8  Y 1  E  (  C  Ofl ,  1  7 1 ,  I  F  A  T  d  ),5,2) 

84 

call 

UYTt  ICON, 175,1  FAT C2> ,5,2) 

85 

call 

BYTE (COM, 179, I  FAT (3), 5, 2) 

86 

CALL 

BYTE  (COM, 1 69, ICAS ( 1 ) ,5,2) 

87 

CALL 

BYTE<CON,173,ICAS(2),5,2) 

88 

CALL 

i>YTE(COr.,177,lCAS<3>,5,2) 

89 

CALL 

BYTE (COM, 184, I  CAT, 1,5) 

90 

CALL 

BYTE ( COM, 043, N AXU06, 1,6) 

91 

call 

BYTE  (  CoM,037,  IiXiDOS,  1  ,6) 

92 

CALL 

aYTt (CON, 193, ISHLVd ) ,4,3) 

93 

CALL 

BYTE  (CO.I,  196,ISHLV(2),4,3) 

94 

CALL 

BYTE  (CUM, 199,  ISHLV (S> ,4,2) 

95 

CALL 

BYTE (CGM, 202,ISHLV(4),4,3) 

96 

CALL 

BYTE(COK,205,ISHLV(5>,4,3) 

97 

CALL 

bYTE (COM, 208, I  SHLV (6) ,4,3) 

98 

CALL 

BYTE(COf-:,223,ISHLV(7),4,3) 

99 

CALL 

BYTE (COM, 191, ISHLP(1),S#2) 

100 

CALL 

dYTE(COI’i,211,ISHLP(8),5,2) 

101 

CALL 

8YT6(C0M,213,IsHLP(3),5,2> 

102 

CALL 

BYTE (COM, 2 15,1 SMLP (4) ,5,2) 

103 

CALL 

bYTL (CON, 2 17, ISHLP (5), 5, 2) 

104 

CALL 

bYTE (COM, 219,ISHLP(6),5,2) 

1  u5 

CALL 

bYTE (COM, 221, ISHLP (7) ,5,2 ) 

106 

IF<  ICAT.6O."75  099".OR.  ICAT.EO." 

107 

ISHLV(1 )=8C0 

108 

I  SHLV ( 2 ) *  250 

109 

1  SHLV  <  3 ) • 1 5 J 

110 

1SHLV(4)*10U 

111 

I SHLV ( 5  >  »40 

112 

ISHLV(6)*20 

75C990" )GQ  To  251 


Table  0-1.  The  modified  SIOAC  post-processor  (continued) 


l SHLV ( 7 ) a  1 Q 
ISHLPC1  )  =  14 
ISHLP(2>*7 
ISHLP<3)*7 
I SHLP ( 4 ) a  1  4 

1  ShlP (51*28 
ISHLP(6)*1  7 
I SHLP ( 7 ) a  1 3 
60  TO  252 

251  I SHLV ( 1 >  *40 
ISHLV<2>«35 
l  SHLV ( 3  >  =30 
I SHLV (41*25 
I SHLV (51=20 
I S  HLV (61*15 
I SHLV (  7 1  *  1  U 
ISHLP(11*2 
ISHLP(2l*5 
ISHLP(31*15 
ISHLP(41*8 

I SHLP ( 5 ) *22 
I  SHLP(61=6 
ISHlP(71=42 

252  CONTINUE 

If  (ICTY.£0."US"1fa0  TO  1017 
00  1011,  JJa1,5Q 

I  f  UCTY  (  J  J  1  .to  ."  "JOO  TO  1017 

1014  I  f (  ICT Y  .EQ.KCT Y < J J ) ) GO  TO  1015 

1015  00  1016  KK*  1 ?7 
I$HLV(MO»XSHLV(X.Kl 

1016  ISHLP(KK1*KSHLP((JJ-11  +  8  +  KK.1 

1017  CONTINUE 

I  f  (  I  FLAG  .66.1  )  60  TO  5  2 
SVCTY»ICTY 
SVRG=IRG 
$ VSUB*  I  SUB 
52  CONTINUE 

CAP*FL0AT(ICAP1 

IF(ICAT.E<3. "75100". OR. ICAT.EQ.” 751000 "1CAP*CAP+1 000 . 
If(ICAT.E6, "75099 ".OR. ICAT.EQ." 750990 ")CAP*CAP«1Q0. 
23  CONTINUE 
I  N  T  =  1 

I  N  F  =  7 

IF  (IFAC.E6.51  I!lT=o 
IF  ClPAC.E«.5l  it.  f  =  1  u 
0  0  191  I  i  =  1  ,  3 
C  A  S  ( 1  3  1  =  J  . 

F AT ( l Jl=u. 

191  CONTINUE 

II  =  1 

44  IFUI.EU.o)  GO  Tu  2C1 

1 F ( ICC ( I  I ) .t0.  ICTY .OK. ICC C I  I 1 .60.  ’ X*’  )  00  To  Ho 

11=11+1 

uO  T  u  UU 

PklwT 

202  font'AT  ( 1  5x  ,//////,•  C  Aliu  XX  MlSilNo'1 
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Table  D-l.  The  modified  SIDAC  post-processor  (continued) 


169 

170 
1  71 

172 

173 
1  74 

175 

176 

177 
176 

179 

180 
181 
162 
163 

184 

185 

186 

187 

188 
169 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 
2 'JO 
2C1 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 
21  7 
218 

219 

220 
221 
222 

223 

224 


CALL  EXIT 

146  00  143  LL*1,7 

143  SHLPCT (LL)*FLOAT( ISHLP(LL) ) /I  00. 

C  ***  THE  CALCULATION  USES  ONLY  THE  PCT  DIST. IN  THE  DATA  H  EC  0  R  0.  hO  «  E  9  c  k 
C  •»*  ISHPCT  HAS  dEEN  READ  IN  ANO  PGM  CAN  d£  MODIFIED 
C  *»»  IT  IS  ASSUMED  THAT  THERE  WILL  6t  ONLY  ONE  PCT  DIST.  FOR  DOTH 
C  *•*  (POSSIBLE)  SETS  OF  PROTECTION  F A C T 0 RS ( I  PR  OF >  . 

C  »**  THE  FORM  OF  THE  CALCULATION  ASSUMED  FoR  THE  SECOND  TIME  PERIOD 
C  ***  IS  1ST  NULT»(1./1ST  PRQT.  FACT.)  ♦  2nD  MULT*(1./2N0  P.F.). 

C  ***  ALSO  IT  IS  ASSUMED  THAT  1ST  PERIOD  UStS  THE  PROT.  FACT.  IN  THE 
C  *»»  DATA  EASE  RECORD, 
c  *  *  * 

C  **»  IPCTF  IS  THE  FLAG  To  USE  NEb  PCT  DIST.  OR  NOT. 

IF(IPCTF(H).E(..0)  GO  TO  144 
DO  145  LL«1,7 

145  SHLPCT(LL)=FLOAT( ISHPCT ( I  I ,LL  > ) /I  00. 

144  IFdI.LT. 6)  GO  TO  147 
DO  148  J  J  =1 , 7 

148  SHLVAL(JJ)»FL0AT(ISHLV(JJ))/1U0. 

60  TO  991 

147  DO  149  K  =  1,7 
PRl=FLOAT (I  PROF  <1 , II/K)) 

PR2*FL0AT(IPROF(2,II,K)) 

IF( IPSOF (1/II.K) .LE.O)  F  R1=FL0AT( ISHLV(K)  > 
IF(IPROF(2#II/K).LE.U)PK2*1 .5  3640 

C  *•*  IF  THE  PROTECTION  FACTOR  IN  THE  DATA  base  RECORD  AND  ThEPROT. 

C  ***  FACT.  IN  THE  CHANGE  CaRu  ARE  bOTM  0/  IT  IS  ASSUiiEO  THAT  THE 
C  •**  PROT.  FACT.  *  1./.65 

XF)ULT(  1/1  l)»FLOAT(  I  MULT  (1  ,  I 1)  )  /10U. 

XMULTC2, I I)=FLOAT< IMULT(2/I I ) ) / 1 00 . 

I  F ( I MU  LT(1,II)  . EU.C. A NO. I  MULT (2,1  I)  .EQ.O)  XMULT(1,II>=1.0 
IF((ICTY.E0."JA”.0R.ICTY.EG."TW”).AND.K.EU.1)HR1*1.4 
lF(ICTY.tC."TH',.AND.2.EG.1)SHLPCT(1 >».3 
IF(ICTY.EG."JA'*.A.ND.K.E«.1)SHLPCT(1>«.4 
I T  T*-24 

N  1  *  ( ( 1 .0  -dd)/ ( AA*24. ) )  +  2 

00  153  NA=1,N1 
00  154  11=1,3 
FATX(I1,NA)*0. 

154  CASX(I1,NA)*G. 

ITT»I TT  +  24 

IF(ITT.Eu.U)»0  TO  152 
I F( ITT .GE .504) GO  TO  151 
ORfc (1 )  =  1  —  ITT** (-.2) 

DRM(2)»ITT**(-.2)-jT f  **(-.2) 

GO  TO  150 

151  uRM(1)=1.-ITT»*(-,2) 

DKK(2)=0. 

GO  TO  150 

152  DRM ( 1 ) =1  . 

DRM(2)=-jTr»*(-.2) 

150  CONTINUE 

FAC1 »DRM( 1 ) »( 1 . /PR1  > 

F  A  C  2  =  X  MUL  T ( 1 ,1  I ) • ( 1  ./PK1  > 

FAC3*XMULT(2,I I )« (1 ./PR  2) 

FAC4*0RM( 2)  • ( F AC2+FAC3) 


83 


Table  D-l.  The  modified  SIDAC  post-processor  (continued) 


225 
22  6 

227 

228 
229 
250 

231 

232 

233 

234 

235 

236 

237 
236 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 
267 
Zba 

269 

270 

271 

272 

273 

274 

275 

276 

277 
27  b 

279 

280 


FAC5*0fiMC1  >+0RM2> 
f AC0*FAC1 +  F AC4 
SHLVALCK)*FAC6/FAC5 
SHVALXCK,NA)»SrtLVALU> 

9563  FORMAT ( 1 X ,1 1 U,6F 1 0.4 ) 

153  CONTINUE 
149  CONTINUE 
991  CONTINUE 
INT»1 
INF*7 

90  999  11*1,3 
DO  999  I2*INT,INF 
DC  999  NA*1,N1 
SMLO (2) *1000. 

S I GL ( 2 ) *200. 

STOR*MAXOOS*SHVAlXCI2,NA) 

IFCSTOR.LE.400.)SMLOC2>=54G. 

IFCST0R.L£.4J0.)SIGLC2)*162. 

IFCSTOR.GT  .40Q..AN0.ST0R.CE.1  300.) SML0C2)*72U. 
1FCSTOH.GT.400..AN9.STOR.CE.1 300. ) SIGLC2>*370. 

FATXC 1 1 ,NA) *FATX(I1,nA)*CUMN(  CST0R-$MLDCI1))/SIoLCI1>>*SHLPCTCI2) 
CASX  C II  ,NA)*CASX ( II ,NA)*CUHNC CSTOR-SMCDC II ) I/SIGCC 1 1 ) >  *SHCPCTCI2 ) 
IFCSTOR.GT. 1000.) J JJCN*j J JCN+1 
IFCSTOR.NE.OJGO  TO  999 
FATXC II, NA)*0. 

CASXCI1,NA)*0. 

999  CONTINUE 

IIN1*  N1-2 
DO  157  11*1,3 

C AS C 1 1 ) *83* C A SX C 1 1 , 1 ) ♦ C 1  .0  - ( 1 1 N1 * AA *24 . )  -GO  )  *  C  A  S  X  C  I  1  ,  N 1  ) 

F A T C I  1 ) *08* F A T X C I  1 , 1 ) ♦ C 1  . 0  - < I  I N 1  *  A A *  24 . )  -8d  )  » F A T X ( I  1  , N 1  ) 

TEMPI *U. 

TEMP2*0. 

00  156  Nb*  2 ,N 1  - 1 
TEMPI • TEMPI *f ATX < I 1,N8) 

156  TEMP2*TEMP2+CASX(I1,N8) 

FATCI1)*FATCI1J*  <  2  4. » A A* TEMPI ) 

1 57  CASCI1 )*CAS(I1 )*C24.*AA«TEMP2) 

00  22  1*1,3 

A  TE  M ( 6 , I) «C AP 

ATEM(5,I)*C1.|J-SPR0BS)*CAP 

atemci,i)*sprobs»cap*fat(i) 

AT£tfC2/l)*(SPR98S  —  F.PROuS) »(1 .U-FAT (  l) ) *CAP 

AT  EMC  3,  I )  *;iHkJ  jS  *  C1.C-FATCI))*CASCI)*CAP 

ATEMC4,  I)  *C  AP-  CATEm  C5,I)*AT£f'1C1,I)*ATEMC2,l)*ATEF’C3,I>> 

IF  CATEMC4,I).LT.O.)  ATEI\C4,I)*0. 

22  Continue 

IF  CIFLAG.Eu.O)  GO  To  4 

IF  C  I  C  T  V  .NE  .  SVC  T  Y  .UK  .  I  hu  .  Nt  .3  VRG  .  Ci«  .  I  SUo  .ME  .  SVSUd)  11*1 
IF  CItTY.Wk.SVCTY.CR.IRu.Me.SVhG.UR.ISUO.Mt.iVSUd)  Go  TO  7 

4  I FL  Ao  *  1 

Oo  5  J 1  *  1 , j 
00  5  J  2  =  1  ,a 
00  5  j  3  *  1  ,3 

AC  I  -i  C  J  1  ,  j  2  ,  J  i  )  *  AL  I  N  C  J  1  ,  J  2  ,  j  5  )  ♦  A  T  t  M  C  J  2  ,  J  3  ) 

5  CONTINUE 
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Table  D-l.  The  modified  SIDAC  post-processor  (continued) 


281 

IF  (IE0F.E</.1)  iO  TO  7 

2  62 

GO  TO  6 

2  83 

7 

ILI..-ILI.J  +  4 

284 

IF  (ILIN.LE.55)  GO  TC  t 

285 

33 

PRINT  69 

286 

69 

FORrIAT  ( 1  H 1  ) 

287 

PRINT  61 

268 

PRINT  02/TITLE1 

289 

PRINT  65*  T I TLE2  /  J MAW 

290 

PRINT  64/  TITLki 

291 

PRINT  05/TITLE4 

292 

PRINT  76 

293 

PRINT  68 

294 

PRINT  66 

295 

61 

FORMAT  ( T2 / *  * »  RED  ON  B LU E 1 /4 0 X/ • PC  PUL A T I ON  ASSESSMENT*) 

296 

62 

FORMAT  ( T2/ *  **  CASE/SCENARIO:  */A20) 

297 

63 

FORMAT  (T2,***  SPEC  INSTR;  * , A20 , 26 X, A 5 ) 

298 

64 

FORMAT  (T2,*»*  SPEC  INSTR:  * / A20 / 2o X, 5 ( 1 H- ) ) 

299 

65 

FORMAT  (T2/'**'/1X/A3/1X/* ASSESSMENT*) 

300 

68 

FORMAT  ( T37, *  FATAL  I  TIES*, S6X,* INJURIES') 

301 

66 

FORMAT!  T22,43<lH-)/2X/4  3C1H-) ,6X, *  RESIDUAL  * / 3X / *  AS S E S S * // / 

302 

S  12X#*  TOT  POP*, 6X,2(*PR0MPT*,8X, 'FALLOUT*, 1UX, ’TOTAL*, 9X>, 

303 

8  2X,*P0P*,7X,*TIME*,/,12X,7(1H-) ,2X,7(1 SC1H-) ,2X),6(1H->,/ / 

304 

8  1X,*REG*,2X,*CTRY*,4X,*MIL*,7X,0(*MIL*,3X, *PCT*,6X>, 

30  S 

8  'MIL*/3X,*PCT' ,4X/*  DATS*///1X/i(1H-)/2X/4(1 H- )  ,  3'X,5  (  1  H - > ,  3X , 

306 

8  7(8(lH-),1X,4(1H-),2X),6(lH->,//> 

307 

I L I N-1 7 

308 

8 

DO  13  *2-1,6 

309 

DO  13  K3» 1 / 3 

310 

ALIN(I1,K2,K3)»ALIN<I1,K2,K3>/MIl 

311 

13 

CONTINUE 

312 

GO  TO  <40, 41, 42), 11 

313 

40 

ENCODE (TOPO S/50)  SVRG/SPAC1/SVCTY 

314 

50 

FORMAT  (A1/A4/A2/*  *) 

315 

GO  TO  75 

316 

41 

ENCODE  (TOPOS/51)  SVRG/SPAC2 

317 

51 

FORMAT  ( A  1 / A  3/ ' ALL  *) 

318 

GO  TO  75 

319 

42 

ENCODE  (TOPOS/57) 

320 

57 

FORMAT  < *WW',2X, 'ALL*  ) 

321 

75 

00  9  13-1/3 

322 

TFA»ALIN(I1/5/l3)+ALIN(I1/1,I3) 

323 

TIN-ALIN (11/2, 13)+ALIN (11,3/13) 

324 

DO  10  K-1/5 

325 

PCT(K)»(ALIN(I1,K,13)/ALIN(I1,6,I3))»100. 

326 

10 

CONTINUE 

327 

PCT (o)*(TFA/ALlN< 1 1/6/13) ) *  1 00  . 

328 

PCT (7)*(TIN/ALIN( I 1/6/13)  )  *100. 

329 

PRINT  1/  TOPOS  ,ALIN( I 1,o/I5) ,ALIN( 11 ,5 ,  I3),PCT<5), 

330 

&  ALIN(I1/1/I3)/PCT(1)/TFA/PCT(6)/ALlN(I1/2/13)/PCT(2), 

331 

8  ALlN(11/3/l5)/PCT(3)/TIN/PCT(7)/ALlNU1/4,13)/PCT(4),ITIME(I3) 

332 

1 

FORMAT  ( 2  X / A3/  1 X / F 8 . 3 / 2 X , 7 ( F 8 . 3 , 1 X , F 4 .  1 , 2 X  )  ,1 X / A  3 ) 

333 

9 

CONTINUE 

334 

IF  (11. EG. 2)  PRINT  76 

335 

IF  (11. EG.  2)  ILIN-ILI N+  1 

336 

PRINT  76 

f 


Table  D-l.  The  modified  SIDAC  post-processor  (continued). 


337 

76 

FORMAT  (/) 

338 

60  11  12=1/6 

339 

DO  11  13*1/3 

340 

ALINCI1/I2/I3) =0 . 0 

341 

1 1 

CONTINUE 

342 

IF  ( 1 1  .  fc  J  .  1  )  SVCTY*ICTY 

343 

IF  (I1.EQ.3.AND.IE0F.EU.1)  60  TO  99 

344 

IF  (11.  £0.2)  SVR6=IRG 

345 

IF  (11. Eu. 2)  60  TO  12 

346 

IF  (11  .60.3)  SVSUb=ISUb 

347 

IF  (I1.E0.3)  I  L  X  fi  =  6  U 

343 

IF  (I1.E0.3)  JNAM  s'RUHAU* 

349 

IF  (11  .to. 3)  60  TO  4 

350 

IF  ( IRG.NE.SVR6.0R.  ISUb.NE.SVSUd.OR  .1  EOF  .Eli  .  1  ) 

351 

IF  (IRG.NE.SVR6.0R  .  I  SUa  .f.E  .  3  VSU6  .  OR  .  I  EOF.EO.  1  ) 

352 

6 u  TO  4 

353 

12 

IF  (ISUb.NE.S6SJb.0R.IEof .£0.1)  11=3 

354 

IF  (ISUb.NE.SUSUb.UR.IECF.EO.I )  60  To  7 

355 

60  TO  4 

356 

96 

I  EOF  =  1 

357 

11*1 

358 

60  TO  7 

359 

99 

STOP 

3oU 

ENO 

86 


APPENDIX  E 

S1DAC  POST-PROCESSOR  RESULTS 


Tables  E-l  to  E-9  present  the  complete  results  of  casualty 
calculations  for  nine  different  stay  time  assumptions.  The  first  seven 
are  traditional  stay  times,  ranging  from  3/4  day  fully  sheltered,  20  1/4 
day  partially  sheltered,  to  21  days  fully  sheltered.  The  last  two  are  the 
Best  and  Worst  Cases  from  the  Attrition  Rate  Model. 

The  column  headings  one  for  the  most  part  self-explanatory.  All 
population  figures  are  in  millions.  The  three  assessment  times  per  table 
are  7  days,  30  days,  and  180  days,  as  indicated  in  the  last  column.  For 
each  table,  the  upper  chart  refers  to  the  urban  population;  the  lower,  to 
the  rural  population. 


Table  E-l.  SIDAC  post-processor  results  -  3/4  day  sheltered,  20  1/4  days  2/3  sheltered. 
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Table  £-5.  SIOAC  post-processor  results  -  7  days  sheltered,  14  days  2/3  sheltered. 
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Table  E-7.  SIDAC  post-processor  results  -  21  days  sheltered. 
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